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ABSTRACT 

We present 2603 spectra of 462 nearby Type la supernovae (SN la), including 2065 previously- 
unpublished spectra, obtained during 1993-2008 through the Center for Astrophysics Supernova Pro- 
gram. There are on average 8 spectra for each of the 313 SN la with at least two spectra. Most of the 
spectra were obtained with the FAST spectrograph at the FLWO 1.5 m telescope and reduced in a 
consistent manner, making this data set well suited for studies of SN la spectroscopic diversity. Using 
additional data from the literature, we study the spectroscopic and photometric properties of SN la 
as a function of spectroscopic class using the classification schemes of Branch et al. and Wang et al. 
The width-luminosity relation appears to be steeper for SN la with broader lines, although the result 
is not statistically significant with the present sample. Based on the evolution of the characteristic 
Si II A6355 line, we propose improved methods for measuring velocity gradients, revealing a larger 
range than previously suspected, from ^ to ~ 400kms~ day" considering the instantaneous ve- 
locity decline rate at maximum light. We find a weaker and less significant correlation between Si li 
velocity and intrinsic B — V color at maximum light than reported by Foley et al., owing to a more 
comprehensive treatment of uncertainties and host galaxy dust. We study the extent of nuclear burn- 
ing and the presence of unburnt carbon in the outermost layers of the ejecta, and report new detections 
of C II A6580 in 23 early-time SN la spectra. The frequency of C ii detections is not higher in SN la 
with bluer colors or narrower light curves, in conflict with the recent results of Thomas et al. Based 
on nebular spectra of 27 SN la, we find no relation between the FWHM of the iron emission feature at 
'^ 4700 A and Ami5{B) after removing the two low-luminosity SN 1986G and SN 1991bg, suggesting 
that the peak luminosity is not strongly dependent on the kinetic energy of the explosion for most 
SN la. Finally, we confirm the correlation of velocity shifts in some nebular lines with the intrinsic 
B — V color of SN la at maximum light, although several outliers suggest a possible non-monotonic 
behavior for the largest blueshifts. 
Subject headings: supernovae: general — surveys 
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Type la supernovae (SN la) have received much atten- 
tion due to their use as distance indicators on cosmolog- 
ical scales and the subsequent discovery of cosmic accel- 
eration (Riess et al. 1998; Perlmutter et al. 1999). Their 
astrophysical nature, however, remains perplexing to this 
day. The standard model of SN la involves the thermonu- 
clear disruption of a carbon-oxygen (C/0) white dwarf 
(WD) star (Hoyle & Fowler 1960; this has been confirmed 
by recent observations of the nearby Type la SN 20 life 
in MlOl; Nugent et al. 2011; Bloom et al. 2012) as it 
approaches the Chandrasekhar mass (A^ch ~ 1.4 M©). 
To do so the WD either accretes material from a non- 
degenerate binary companion (the "single-degenerate" 
scenario), or merges with another WD star (the "double- 
degenerate" scenario; Iben & Tutukov 1984; Webbink 
1984). Both scenarios constitute plausible progenitor 
channels, and can accomodate variations in the total 
mass of the binary system (i.e., sub- and super- Mch! see 
Howell 2011 for a recent review). 

The favored explosion mechanism involves the tran- 
sition from a turbulent subsonic burning front propa- 
gating via thermal conductivity (known as a defiagra- 
tion; e.g., Nomoto et al. 1984) to a supersonic detona- 
tion propagating via a strong shock (Khokhlov 1991). 
The defiagration phase synthesizes iron-peak elements, 
of which ^^Ni is needed to power the light curve through 
the ^^Ni-^^^Co— !>'^^Fe decay chain. It also pre-expands 
the WD so that the subsequent detonation burns the left- 
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over C/0 fuel at lower densities and synthesizes sufficient 
amounts of intermediate-mass elements (IME; e.g., Si, S, 
Ca) at high expansion velocities, needed to reproduce ob- 
served spectra (e.g.. Branch et al. 1982). Recent multi- 
dimensional simulations highlight the importance of hy- 
drodynamical instabilities during the deflagration phase 
(e.g., Gamezo et al. 2005; Ropke & Niemeyer 2007), 
and provide a physical basis for the deflagration-to- 
detonation transition (DDT; e.g., Woosley et al. 2009), 
although this DDT is artificially triggered based on pre- 
defined criteria (see, e.g., Kasen et al. 2009; but see 
Poludnenko et al. 2011 for 3D simulations of spontaneous 
DDTs in unconfined media). 

Since the first spectrophotomctric measurements of a 
SN la (SN 1972E; Kirshner et al. 1973), detailed ob- 
servational studies of individual objects have revealed 
a wide range of photometric and spectroscopic proper- 
ties, including luminous 1991T/1999aa-like events (e.g., 
Filippenko et al. 1992b; Garavini et al. 2004), faint 
1991bg-like SN (e.g., Filippenko et al. 1992a; Leibundgut 
et al. 1993; Garnavich et al. 2004; Taubenberger et al. 
2008), peculiar 2002cx-like events (e.g., Li et al. 2003; 
Phillips et al. 2007), and high-luminosity SN la specu- 
lated to originate from super-Chandrasekhar-mass pro- 
genitors (e.g., HoweU et al. 2006; Scalzo et al. 2010). This 
observed diversity contrasts with the apparent homo- 
geneity of the SN la class, the majority of which seems to 
obey an empirical relation between the peak luminosity 
and the width of the light curve (the so-called width- 
luminosity relation, or WLR; Pskovskii 1977; Phillips 
1993), needed to normalize SN la luminosities for pre- 
cise distance measurements. Several large sets of optical 
SN la light curves have already been published (Hamuy 
et al. 1996; Riess et al. 1999; Jha et al. 2006; Hicken et al. 
2009; Contreras et al. 2010; Ganeshalingam et al. 2010; 
Stritzinger et al. 2011), but until the recent publication 
of 432 spectra of 32 SN la by Matheson et al. (2008), 
there had been no large homogeneous spectroscopic data 
sets of nearby SN la available. 

The supernova group at the Harvard-Smithsonian Cen- 
ter for Astrophysics (CfA) initiated an ambitious obser- 
vational program in 1993 (PI: R. P. Kirshner) using the 
telescopes at the Fred Lawrence Whipple Observatory 
(FLWO) on Mount Hopkins, Arizona. The aims of the 
CfA SN Program are to obtain well-sampled light curves 
and multi-epoch spectroscopy of nearby {z < 0.015) su- 
pernovae of all types, to study their astrophysical nature 
as well as provide a homogeneous low-redshift sample for 
cosmological analyses using SN la. We follow supernovae 
announced in lAU circulars^^, discovered for the most 
part by the Lick Observatory Supernova Search (LOSS) 
using the robotic Katzman Automatic Imaging Telescope 
(KAIT; Filippenko et al. 2001), and by dedicated ama- 
teur astronomers. In return we provide spectroscopic 
classifications for a large fraction of supernova discover- 
ies: during 1997-2007, we classified over 40% of all nearby 
SN visible from Mount Hopkins and announced our re- 
sults in 499 lAU circulars. These prompt classifications 
are important to the community, and serve as the basis 

^^ Starting in 2008 wc also followed a few supernovae dis- 
covered by the SNfactory collaboration (Aldering et al. 2002) 
and reported on the SNfactory Supernova Discoveries webpage, 
http://snfactory.lbl.gov/snf/opcn_access/snlist.php. 



for triggers of other supernova follow-up programs (e.g., 
by the Swift SN group; Milne et al. 2010). 

We were not able to follow up all SN announced in cir- 
culars, and so had to prioritize our follow-up strategy. In 
a nutshell, we assigned the highest priority to (1) SN la 
discovered before maximum light (we sometimes relaxed 
this requirement to follow up peculiar objects) (2) bright 
Type Ib/c supernovae, to study both individual objects 
(e.g., SN 2006aj, Modjaz et al. 2006; SN 2008D, Mod- 
jaz et al. 2009) and the connection of broad-line SN Ic 
with gamma-ray bursts (e.g., Modjaz et al. 2008), and 
(3) bright Type IIP supernovae discovered within '^ 2 
weeks from explosion for distance determinations via a 
variant of the Expanding Photosphere Method (Kirshner 
& Kwan 1974) [e.g., Dessart et al. 2008]. 

We have already published several large collections 
of optical SN la light curves obtained with the FLWO 
1.2m telescope (CfAl: 22 BVRI light curves, Riess 
et al. 1999; CfA2: 44 UBVRI light curves, Jha et al. 
2006; CfA3: 185 UBVRr'Ii' light curves, Hicken et al. 
2009), as well as a sample of 21 near-infrared (NIR; 
JHKs) light curves obtained with the robotic PAIRI- 
TEL 1.3 m telescope (Wood-Vasey et al. 2008). Spectra 
of individual SN la have been published in several papers 
(SN 1998aq, Branch et al. 2003; SN 1998bu, Jha et al. 
1999b; SN 1999by, Garnavich et al. 2004; SN 2001ay, 
Krisciunas et al. 2011a; SN 2001eh and SN 2001ep, Sauer 
et al. 2008; SN 2002cx, Li et al. 2003; SN 2005cf, Wang 
et al. 2009b; SN 2006bt, Foley et al. 2010; SN 2006gz, 
Hicken et al. 2007a), including the large sample of 432 
spectra of 32 SN la (including SN 1998aq and SN 1998bu) 
published by Matheson et al. (2008). A study of the use 
of spectroscopic indicators in determining distances to 
SN la in the Hubble flow was published by Blondin et al. 
(2011b), based on 47 spectra of 26 SN la (including 15 
spectra of 9 SN la from the Matheson et al. 2008 sample). 

Here we present the second large release of SN la spec- 
tra from the CfA SN Program taken between October 
1993 and August 2008. This sample consists of 2603 
spectra of 462 SN la, of which 2065 spectra of 406 SN la 
are published here for the first time. It includes new re- 
ductions of spectra taken during 1994-1997 of SN la that 
were part of the CfAl light-curve sample. Unlike Mathe- 
son et al. (2008), we include spectra of all SN la^^ taken 
during this 15-year period, regardless of sampling and 
whether or not a well-calibrated light curve was available. 
The main reason for doing so is that a large fraction of 
the 149 SN la for which we only have a single spectrum 
were classified by us, and we wish their spectra to be 
publicly available. Another reason is that other groups 
might have well-sampled light curves of objects for which 
we were unable to determine a time of maximum light, 
adding value to the "phaseless" ^^ spectra presented here. 

The resulting sample spans a large range of SN la 
properties and is ideal for studies of spectroscopic di- 
versity, which is the aim of this paper. In § 2 we present 
the spectroscopic observations and data reduction tech- 
niques. Section 3 comments on the general properties of 
the CfA sample of 462 SN la. In § 4 we present the spec- 

^* except for 11 spectra of 10 SN la with too low signal- to- noise 
ratio (S/N) and excessive host-galaxy contamination, which were 
not used for spectroscopic classifications. 

^^ Throughout this paper the spectral phase corresponds to the 
number of rest-frame days from S-band maximum light. 
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troscopic elassification schemes of Branch et ah (2006) 
and Wang et ah (2009a), and study the properties of 
SN la in these different subclasses. Section 5 focuses 
on the characteristic Si ii A6355 Une of SN la, studying 
its velocity evolution and relation to intrinsic color. In 
§ 6 we analyze the extent of nuclear burning and the 
presence of high- velocity features and unburnt carbon in 
early-time spectra. Last, § 7 examines the relation be- 
tween the FWHM of nebular iron lines and luminosity, 
and between nebular line shifts and intrinsic properties 
of SN la at maximum light. Conclusions follow in § 8. 

All 2603 spectra presented in this paper will be made 
publicly available through the CfA Supernova Archive^^. 

2. OBSERVATIONS AND DATA REDUCTION 

In this section we summarize our observational strat- 
egy and data reduction methods. These are essentially 
identical to the ones presented by Matheson et al. (2008), 
but we highlight a few differences and provide additional 
information when needed. 

2.1. Observations 

As with the sample of 432 spectra published by Math- 
eson et al. (2008), the majority of spectra presented here 
(2447 out of 2603 spectra, or ^ 94%) were obtained with 
the FAST spectrograph (Fabricant et al. 1998) mounted 
on the Tillinghast 1.5 m telescope at FLWO. The FAST 
spectrograph has been operational since January 1994 
and the first SN la spectrum was one of SN 1994D taken 
on 1994 March 10 (all dates are given in UT). The obser- 
vations were carried out in queue-scheduled mode, for the 
most part by two professional observers (P. Berlind and 
M. L. Calkins, who observed '- 72% of the FAST spectra 
presented here), as well as by CfA personnel. A total of 
79 individual observers contributed to the FAST SN la 
sample presented here. During 1997-2008, typically 2- 
3 spectra were taken each night FAST was scheduled on 
the FLWO 1.5 m, namely ^ 20 nights per month, exclud- 
ing August which corresponds to the annual shutdown 
at FLWO during the monsoon season. During 1994- 
1996, typically a single spectrum was taken in any given 
night. Observational details of the spectra are given in 
Table Al. 

The usual setup for observations with FAST consisted 
of a ruled grating with 300 lines per mm and a 3" slit, 
yielding a typical FWHM resolution of 6-7 A over a wave- 
length range of ^ 3700 to ^ 7500 A. From September 
2004 onwards (starting with observations of SN 2004dt), 
we changed the standard setup, extending it down to 
^3500 A to cover the entire Ca iiH&K absorption pro- 
file. For a few SN la (SN 1996ai, SN 1999by, SN 1999dq, 
SN 2001ch, SN 2001ep, and SN 2002bo) we requested 
additional observations with a different grating tilt to 
extend the wavelength range beyond 9000 A. Other pro- 
grams in the FAST queue would sometimes request dif- 
ferent instrument setups, either with narrower slits (1.5" 
and 2") or with a higher-resolution grating (cf. spectra 
of SN 1995ac, SN 1998aq, SN 1998bu, SN 1999by, and 
SN 1999cl). For most of the observations obtained dur- 
ing 1994-1998, the slit was oriented at a position angle 
of 90°. From late 1998 onwards, the slit was generally 
positioned at the parallactic angle (unless the object was 

^^ http;//www.cfa. harvard. cdu/supcrnova/SNarchivc. html 



at airmass < 1.1) so as to minimize the effects of atmo- 
spheric dispersion (Filippenko 1982). Table Al gives the 
slit position angle for each spectrum as well as the abso- 
lute difference, |A$|, with the actual parallactic angle. 

Additional spectra were obtained during classically- 
scheduled nights at the MMT Observatory 6.5 m tele- 
scope with the Blue Channel (131 spectra) and Red 
Channel (3 spectra) spectrographs (Schmidt et al. 1989). 
A number of different spectrograph settings were used, 
yielding FWHM resolutions ranging between ^ 3 and 
^ 13 A, with a wavelength range typically extending be- 
low 3500 A and beyond 8000 A. A few spectra were taken 
with two settings of the Blue Channel spectrograph with 
non-overlapping wavelength ranges. This only concerns 
the five MMT spectra taken on 1994 June 12 (SN 1994D, 
SN 1994M, SN 1994Q, SN 1994S, and SN 1994T), all of 
which have ^ 400A-wide gaps (~6150-6550A) between 
the two portions of the spectra. These spectra have been 
clearly marked in Table Al. 

Last, a few spectra were obtained with the Magellan 
6.5 m Clay (-hLDSS-2/LDSS-3; 21 spectra) and Baade 
(-hIMACS; one spectrum of SN 2003kf) telescopes. The 
FWHM resolution varies between '^ 9 and ~ 18A for 
LDSS-2, 9-12 A for LDSS-3, and 5-6 A for IMACS. The 
spectra taken with the Magellan 6.5 m telescopes typ- 
ically do not reach bluer wavelengths than the FAST 
spectra, but they generally extend beyond 9000 A. 

2.2. Data Reduction 

The data reduction methods are the same as those pre- 
sented by Matheson et al. (2008), and we refer the reader 
to that paper for complementary information. The FAST 
data were all reduced in the same consistent manner. 
This also applies to FAST data from 1994-1997 that were 
entirely rc-rcduced for the purposes of this paper (this 
enabled the recovery of two spectra of SN 1997do and 
SN 1997dt accidently omitted from the data set pub- 
hshed by Matheson et al. 2008). The rcsuUing FAST 
sample of 2447 spectra makes this by far the largest ho- 
mogeneous SN la spectroscopic dataset to date. 

We used standard routines in IRAF^^ to correct the 
CCD frames for overscan. The bias frames are suffi- 
ciently uniform so that we do not subtract them to avoid 
introducing additional noise. We did not correct for 
dark current as it is typically negligible with FAST. As 
noted by Matheson et al. (2008), however, a few spectra 
are affected by dark-current problems after UV flash- 
ing, resulting in a small emission feature at ^ 7100 A 
(observed) and increased noise in the red portion of 
the spectrum. For a handful of spectra taken in 1998 
(SN 1998CS, SN 1998de, SN 1998ec, and SN 1998eg), the 
S/N degradation is large enough that wc have simply 
trimmed off affected portions of the spectra (these have 
been clearly marked in Table Al). The CCD frames 
are then flat-fielded using a combined normalized flat- 
field image. One-dimensional spectra are then optimally 
extracted using the algorithm of Home (1986) as imple- 
mented in the IRAF apall package. These spectra are 
then wavelength calibrated using HeNeAr lamps taken 

^^ IRAF is distributed by the National Optieal Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement writh 
the National Science Foundation. 
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immediately after (occasionally before) each exposure. 
The same procedure was applied to spectra taken with 
the MMT and Magellan 6.5 m telescopes, although we 
subtracted a combined bias frame before the flat-fielding 
stage. Subsequent reduction steps described below are 
common to all telescope/instrument combinations. 

We use our own set of routines in IDL to flux 
calibrate the extracted one-dimensional wavelength- 
calibrated spectra. Aside from the actual flux calibra- 
tion, these routines also apply small adjustments to the 
wavelength calibration based on night-sky lines in the 
SN frames, and apply a heliocentric correction. The 
spectrophotometric standard stars observed on the same 
night (see Table Al) are also used to remove telluric ab- 
sorption features from the spectra (see, e.g.. Wade & 
Home 1988; Matheson et al. 2000). We did not re-reduce 
the 37 MMT spectra from 1993-1997 which were flux cal- 
ibrated using standard procedures in IRAF. 

Approximately 50% of the spectra listed in Table Al 
were cross-calibrated with two spectrophotometric stan- 
dard stars of different colors to overcome the impact of 
second-order light contamination. Prior to the FAST re- 
furbishment in August 2003, the spectrograph blocked 
blue light such that the FAST spectra until SN 2003gq 
do not suffer from this effect (FAST spectra taken be- 
tween September 1997 and July 2003 were usually cali- 
brated with a single standard star). Starting in Septem- 
ber 2003, two standard stars were systematically used 
for the flux calibration. Our new reductions of the 1994- 
1997 data also uses two standard stars for consistency 
with the 2003-2008 data set, although the impact on the 
resulting flux calibration is minor. 

Figure 1 illustrates our flux calibration technique with 
two standard stars on the FAST spectrum of SN 2008bf 
at 1 d past S-band maximum. On the same night as a 
given SN observation, we observed a relatively blue stan- 
dard star (here Feige 34, sdO spectral type; top panel) 
and a relatively red standard star (here HD 84937, sdF 
spectral type; middle panel). The blue standard yields 
a more accurate flux calibration in the blue due to a 
greater number of counts and the absence of a Balmer 
jump which complicates the calibration around 4000 A, 
while the red standard yields a better calibration in the 
red where it does not suffer as much from second-order 
contamination. One clearly sees that the blue standard 
Feige 34 is poorly calibrated by the red standard blue- 
ward of ^ 4000 A (top panel; red line), while applying the 
blue standard to calibrate the red standard HD 84937 
results in a flux dcflcit redward of ~ 6000 A {middle 
panel; blue line). The SN spectrum is calibrated us- 
ing both standard stars (bottom panel), and the two re- 
sulting spectra arc combined in a ~ 100 A- wide region 
around ^ 4500 A (i.e., slightly redward of the Balmer 
jump). Blueward (redward) of this region, the spectrum 
calibrated using the blue (red) standard is used. The 
impact on the final spectrum (bottom panel; black line) 
may not seem spectacular, but we clearly recover fiux 
in the Ca ii H&K absorption region and redward of the 
Si II A6355 line. 

We also illustrate the removal of atmospheric absorp- 
tion features, here the O2 B-band around 6880 A. Each 
standard star is used to derive a normalized telluric spec- 
trum (^top and middle panels; black line) which is then 
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Fig. 1. — Illustration of our flux-calibration procedure with one 
"blue" (Feige 34; top panel) and one "red" (HD 84937; middle 
panel) spectroscopic standard star, applied to the FAST spectrum 
of SN 2008bf at one day past _B-band maximum {bottom panel). 
In all three panels, the blue (red) spectrum corresponds to a flux 
calibration using the blue (red) standard star. The gray shaded 
area highlights the wavelength interval used to combine the blue- 
and red-calibrated spectra of SN 2008bf. Blueward (redward) of 
this interval, the blue-(red-) calibrated spectrum is used to generate 
the final spectrum {black line). Both standard stars are also used 
to generate normalized telluric absorption spectra {black line in 
top and middle panels) which serve to remove these features (here 
the atmospheric B-band) in the final SN spectrum {see inset in 
bottom panel; black line). 

divided out from the SN spectrum after appropriate scal- 
ing (ratio of airmasses raised to some power a « 0.6; see 
Wade & Home 1988). The red standard star is more 
effective in removing these features, since the blue stan- 
dard suffers from second-order contamination in these at- 
mospheric bands. In this particular example the B-band 
leaves a small imprint on the SN spectrum, but the inset 
in the lower panel of Fig. 1 shows that this technique is 
successful in removing the unwanted absorption. 

We can check the accuracy of the relative flux cali- 
bration of our spectra (no attempt was made to put the 
spectra on an absolute flux scale) by comparing the B — V 
color derived from photometry with that derived directly 
from the spectra (see Matheson et al. 2008, their Fig. 4). 
We interpolated the corresponding B- and F-band light 
curves at the time each spectrum was taken, unless the 
difference was larger than 3 d (the interpolated measure- 
ments were visually cross-checked). FAST spectra taken 
before September 2004 do not extend to the blue edge 
of the B filter (3600 A). We have run simulations based 
on the SN la spectral template of Hsiao et al. (2007) 
that show that even for spectra extending only down to 
3750 A the error on the inferred B — V color is less than 
0.005 mag, so we include ah the 1994-2008 FAST spec- 
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Fig. 2. — Comparison oi B — V color derived from spectra with 
that derived from photometric measurements (see Matheson et al. 
2008, their Fig. 4). We show measurements for spectra with phases 
less than 20 d past B-band maximum ( circles) , for which the dif- 
ference between the position angle of the spectroscopic slit and the 
parallactic angle, noted lA^I, is less than (black) or greater than 
(blue) 10°. Spectra taken at low airmass (sec z < 1.1) are included 
in the |A<I>| < 10° measurements. The stars correspond to spectra 
at phases greater than 20 d past _B-band maximum. In each case we 
indicate the magnitude scatter about the 1:1 relation {dashed line). 
The +63 d spectrum of SN 2005M has low S/N, SN 1996ai is the 
most highly-reddened SN la in our sample {E(B — V) sd 1.8 mag), 
and the spectra of SN 2005cc are contaminated by host-galaxy 
light. 

tra in the comparison. The resulting 1128 B — V mea- 
surements are displayed in Fig. 2. For SN spectra at 
phases less than 20 d past i?-band maximum light taken 
within 10° of the parallactic angle (or at low airmass: 
secz < 1.1; filled circles)^ the scatter around zero differ- 
ence is ^ 0.08 mag, slightly larger than the ~ 0.06 mag 
scatter found by Matheson et al. (2008). The reason for 
this larger scatter is the inclusion of low-S/N spectra in 
our sample (S/N < 10 per pixel). Three significant out- 
liers (~ 1 mag difference m B — V color) correspond to 
spectra of the 2002cx-like SN 2005cc that are contami- 
nated by host-galaxy light. 

The scatter is '^ 2 times larger for spectra not observed 
at the parallactic angle (|A$ > 10° |), regardless of phase 
[open circles and open stars). The scatter is also signifi- 
cantly larger for SN spectra at phases greater than 20 d 
past maximum observed at the parallactic angle {filled 
stars). Possible reasons for this were already noted by 
Matheson et al. (2008): the spectra become dominated 
by prominent emission features, giving rise to systematic 
errors when multiplied by a filter that is not precisely 
matched to the photometry. The spectra are also fainter, 
increasing the impact of host-galaxy contamination. 

We show example spectral series in Fig. 3. Plots of all 
SN la from the CfA SN Program will be made available 
alongside the actual data on the CfA SN Archive website. 

3. SAMPLE PROPERTIES 

We briefly summarize the properties of the SN la sam- 
ple presented in this paper in terms of phase, redshift, 
and Ar?T,i5(i?) distributions, as well as the typical num- 
ber of spectra per SN. We cross-checked the classifica- 
tions of all supernovae in our sample with those reported 
in circulars using the SNID code of Blondin & Tonry 



(2007). Three SN la have no previously reported classi- 
fications: SN 2001es, SN 2003de, and SN 2004cu. One 
SN fa (SN 1999bv) was mis-classified as a "probable Type 
Ib/c" supernova by Jha et al. (1999a), although its spec- 
trum shows it to be a Type fa SN at ^ 3 months past 
maximum light (as noted by Matheson et al. 2001). All 
other SN in our sample are classified as Type la in cir- 
culars. 

We were able to determine a reliable time of S-band 
maximum (tmax(-S); see Table 1) for 271 out of 462 SN fa, 
based for the most part on photometry from the CfA 
(Riess et al. 1999; Jha et al. 2006; Hicken et al. 2009), 
and from the Lick Observatory Supernova Search (LOSS; 
Ganeshalingam et al. 2010). We rely on a number of 
light-curve fitters for this purpose: MLCS2k2 (Jha et al. 
2007), SALT2 (Guy et al. 2007), SNooPy (Burns et al. 
2011a), and the BayeSN statistical model of Mandel et al. 
(2009, 2011). When more than one estimate for tmax(-B) 
is available, we simply take the mean value. For a few SN 
we used values from the literature, in particular results 
from SNooPy fits to light curves from the Carnegie Su- 
pernova Program (CSP; Contreraset al. 2010; Stritzinger 
et al. 2011) by Burns et al. (2011b) and Stritzinger et al. 
(2011). Estimates for six SN la (SN 2001gc, SN 2002es, 
SN 2007fb, SN 2007fs, SN 2007kk, and SN 2008ae) 
were provided by the Berkeley SN group based on LOSS 
light curves (Ganeshalingam, Li, & Silverman 2011; priv. 
comm.), while imax(-B) values for 4 of the 5 SNF SN la 
in our sample were provided by the SNfactory collabo- 
ration (Bailey 2011; priv. comm.). We thus have 191 
SN la with no reliable time of maximum (24 of which 
have B-band light curves). While this represents a size- 
able fraction of the total number of individual SN la in 
our sample, it is only a modest percentage of the total 
number of spectra (397 in 2603 spectra, or ^ 16%). 

The resulting phase distribution of the 2206 spectra for 
which we have a reliable tmax(S) estimate is shown in 
the left panel of Fig. 4 out to 150 d past maximum. Also 
shown is the corresponding distribution for the sample of 
spectra published by Matheson et al. (2008), including 
9 spectra that were accidently omitted from their sam- 
ple (see Table Al). There are 43 spectra at phases less 
than — lOd, and 381 spectra within 3d from maximum 
light. Many of the SN la in this sample have observations 
starting before maximum light. The right panel of Fig. 4 
shows the distribution of SN phases of the first spectrum. 
There are 23 SN la for which the spectroscopic follow- 
up was initiated earlier than — lOd, and 168 SN la with 
spectra starting prior to maximum light. We also have 
21 spectra of 12 SN la with spectra taken later than 150 d 
past maximum (see Table 2). 

A clearer measure of the scientific potential of this sam- 
ple is the number of SN with at least TV spectra. Figure 5 
shows the corresponding distribution, for all 462 SN la 
in our sample (solid line) and for the 191 with no time of 
maximum ( dotted line) . There are 92 SN la for which we 
have at least 10 spectra, three of which have no reliable 
imax(5) (SN 2002bz, SN 2003fd, and SN 2003iu). There 
are 20 SN la with at least 20 spectra, and 4 objects with 
over 30 spectra (SN 1998bu, SN 2001V, SN 2002bo, and 
SN 2007af). There are on average 8 spectra for each of 
the 313 SN fa with at least two spectra. 

The redshift distribution of our sample is shown in 
Fig. 6. The median redshift is - 0.023, and 361 SN la 
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Fig. 3. — Example SN la spectral series from the CfA SN Program. The flux units are />, (erg s~^ cm~^ A~^) that have been normalized 
and then additive offsets applied for clarity. The zero-flux level for each spectrum is marked with an extension on the red edge. The 
wavelength axis is corrected for the recession velocity of the host galaxy. The number associated with each spectrum indicates the phase 
in (rest-frame) days from B-band maximum. Spectra with low S/N have been binned; they are indicated with an asterisk appended to 
the phase label. Shown are spectra of two "normal" SN la in the top row (SN 2007af and SN 2007co); one 1991bg-like (SN 2005ke), one 
1991T-like (SN 2007S; the -|-59d and +64 d spectra are strongly contaminated by the host-galaxy), and one 2002cx-likc (SN 2008A) SN la 
in the bottom row. 



are at z > 0.015 (nearby Hubble flow). The highest- 
redshift SN la in our sample, SN 1996ab at z w 0.123, 
is part of the CfAl light-curve sample of Riess et al. 
(1999). Nearly all redshifts are from the corresponding 
host galaxy as reported in the NASA/IPAC Extragalac- 
tic Database (NED; ~ 40% are from the Updated Zwicky 
Catalog of Falco et al. 1999). For a few objects there is 
no NED redshift for the host galaxy, and we then de- 
termine our own redshift based on narrow emission or 
absorption lines from the host galaxy (either in the SN 
spectrum itself or on a separate host-galaxy spectrum; 
see Table 3), using the IRAF RVSAO package of Kurtz 
& Mink (1998). When no such lines arc visible in the SN 
spectrum, we resort to cross-correlation with a library of 
supernova spectra using SNID, and report the median 
redshift based on the top-5 matches. For two SN la host 
galaxies the redshift reported by NED is in error: MCG 
+07-36-33, the host of SN 2003fa, has a NED redshift 



of 1800 km s (de Vaucouleurs et al. 1991) where we 
measure z = 0.0404 (cz = 12125 km s~^ sec Table 3) 
based on five emission lines (including [O ii] A3727 and 
Ha). The host of SN 2002es, UGC 2708, has a NED 
redshift of 8523 km s~ (Monnier Ragaigne et al. 2003) 
reported on the main page, where the correct redshift of 
5394 km s~ (de Vaucouleurs et al. 1991) is reported on 
the redshift subpage (this is confirmed by the Berkeley 
SN Group using an optical spectrum of UGC 2708; Li 
2011, priv. comm.). All other NED redshifts were found 
to be consistent with our own estimates. Four of the five 
SN la in our sample discovered by the SNfactory occur in 
anonymous hosts. We derive their heliocentric redshifts 
based on the CMB-frame redshifts given in Table 2 of 
Bailey et al. (2009) and the coordinates as reported on 
the SNfactory Supernova Discoveries webpage. 
The distribution of the light-curve width parameter 
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TABLE 1 
Basic SN parameters 



SN 



4B)^ 



Ami5(Br 



Mb" 



B -V 



_B" 



Phot. Ref.s 



0.0503 


49269.2 (1.2) [Jha07] 


1.19 (0.10) [Phi99] 


0.0190 


h 




0.0029 


49432.5 (0.1) 


B,M,SA,SN 


1.37 (0.03) 


B 


0.0232 


49474.5 (0.9) 


B,M,SA,SN 


1.26 (0.06) 


B 


0.0295 


49496.2 (1.1) 


Jha07] 


1.03 (0.10) 


Phi99] 


0.0151 


49518.3 (0.3) 


B,M,SA,SN] 


0.94 (0.06) 


B| 


0.0346 


49514.4 (0.4) 


M,SN] 


1.51 (0.08) 


SN] 


0.0043 


49685.5 (0.8) 


B,M,SA,SN 


0.96 (0.04) 


B 


0.0067 


49768.7 (0.1) 


B,M,SA,SN 


1.05 (0.03) 


B 


0.0120 


49775.2 (0.5) 


B,M,SA,SN] 


1.16 (0.05) 


B] 



-19.21 (0.08) -0.11 (0.06) 
-19.55 (0.09) --0.12 (0.05) 





CfAl 




CfAl 




Jha07 


0.10 (0.06) 


CfAl 




CfAl 


0.11 (0.05) 


CfAl 




CfAl 




CfAl,Jha07 




CfAl,Jha07 




CfAl 



References. — (BSNIP) Berkeley SN la Program (2011, priv. comm.); (Buf05) Bufano et al. (2005); (Burll) Burns et al. (2011b); 
(CfAl) Riess et al. (1999); (CfA2) Jha et al. (2006); (CfA3) Hicken et al. (2009); (ConlO) Contreras et al. (2010); (GanlO) Ganeshalingam 
et al. (2010); (Ganl2) Ganeshalingam et al. (2012); (Hic07) Hicken et al. (2007a); (Jha07) Jha et al. (2007) (and references therein); (KasOS) 
Kashwal et al. (2008); (KrilO) Krisciunas et al. (2011b); (Phi99) PhilUps et al. (1999); (Phi07) Phillips et al. (2007); (Pri09) J. L. Prieto 
(2009, priv. comm.; based on the method of Prieto et al. (2006); (ScalO) Scalzo et al. (2010); (SNF) SNfactory collaboration (2011, priv. 
comm.); (Strll) Stritzingcr et al. (2011); (TauOS) Taubenberger et al. (2008). 

Note. — Table 1 is published in its entirety in the electronic edition of The Astronomical Journal. A portion is shown here for guidance 
regarding its form and content. 

^Heliocentric redshift as reported in NED (except for the redshifts given in Table 3). 
MJD at B-band maximum light. The la error is given in between parentheses. The fitters used to determine imax(S) are given in between 
square brackets: B^BaycSN (Mandcl et al. 2009); M^MLCS2k2 (Jha et al. 2007); SA^SALT2 (Guy et al. 2007); SN^SNooPy (Burns et al. 
2011a). References to estimates taken from the literature are given below. In the case of multiple imax(^) estimates the error is given as the 
standard deviation of the individual values. For cases where no error information was available we assume a typical 0.5 d uncertainty. 

'^ d^mi^{B) decline rate in magnitudes. The Icr error is given in between parentheses. The fitter used to determine Ami5(B) is given in the 
square brackets (sec above). We assign a minimum 0.05 mag uncertainty to /S.m-i^{B) estimates with SNooPy (which typically yield a formal 
~ 0.01 mag uncertainty). References to estimates taken from the literature are given below. 

Intrinsic absolute peak B-band magnitude from BaycSN (assuming Hq — 72kms~ Mpc~ ) for SN la with zcm.^ > 0.01 and Av < 1 mag. 
^Intrinsic B — V color at B-band maximum light from BaycSN for SN la with ZQM^n ^ 0.01 and Ay < 1 mag. 

Intrinsic B"^^^ — y^^^ pseudo-color from BayeSN for SN la with zcmb ^ 0.01 and Ay < 1 mag. 
^References for published photometric data used to determine iinax(^) s^nd. AmiQ{B) (see below). References to "Jha07'' arc meant to point 
the reader to references given in Table 1 of Jha et al. (2007). 

Jha et al. (2007) give tmax(S) — 49288.08 it 1.19, but the fit is uncertain (first photometric point > +10d past S-band maximum). 
Mha et al. (2007) give imax(-B) ^ 50630.25 ± 0.98 and Ami^{B) ^ 1.02 ±0.10, but the fit is uncertain (first photometric point > +10d past 
S-band maximum), 
jjha et al. (2007) give tmax(-B) ^ 50840.57 ± 2.02, but the fit is qualified by them as "poor". 



100 



50 



30 




' 


1 25 




- 


1 20 

i 15 


_n 


- 


„ 10 


^ 


LP 




^. . n 



60 



80 



100 120 140 



-10 



10 20 30 40 

Rest-frame phase [days] 



50 



60 



25 


- p 


\//A Matheson et at (2008) ^ 


20 


- 


- 


15 
10 


- 




J 




\\ 


r 


-n 


- 


5 



L 1 L 


7n / //r-^/} . • * p-71 ( — in 


\f\^\ 



-10 10 

Phase of first spectrum [days] 



20 



Fig. 4. — Distributions of SN la for which we have a reliable estimate for the time of maximum light. Lejt: Distribution of the number 
of spectra as a function of phase out to 60 d past maximum light. The hatched histogram corresponds to the SN la sample published by 
Mathcson et al. (2008). The inset shows spectra at phases in the range [-1-60,-1-150] d. Right: Number of SN la vs. the phase of the first 
spectrum. 

with previously-published values). For the remaining 39 
objects we take A77ii5(i3) values published in the litera- 
ture (mostly from Jha et al. 2006, Ganeshalingam et al. 
2010, Burns et al. 2011b and Stritzinger et al. 2011). 
Our sample covers a large range in A?7ii5(i3), from lumi- 
nous SN la with Ami5{B) w 0.7 mag to faint ones with 
/S.mi^{B) > 2.0 mag. There are 26 SN in our sample with 



ATOi5(i?)^* is shown in Fig. 7 (see also Table 1). Of 
the 271 SN la in our sample with imax(-B) information, 
247 have a reliable A77ii5(i3) estimate. For most SN we 
use the ATOi5(i?) inferred from our own BayeSN (185) 
or SNooPy (23) fits (we checked these were consistent 

^* The difference in B-band magnitude between maximum light 
and 15 d after maximum (Phillips 1993). 
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TABLE 2 

CfA SN Ia spectra at 150 d past 
maximum or later (see table al 
for more observational details). 



TABLE 3 

SN REDSHIFTS NOT TAKEN FROM NED. 



SN 


Phase (d) 


1994D 


+611 


1994ae 


+152,4-366 


1995D 


+275, +283 


1998aq 


+209, +229, +239 


1998bu 


+178,+189,+207,+216,+242 


2000CX 


+154,+181 


2002bo 


+310 


2002CX 


+308 (low S/N) 


2003kf 


+397 


2006X 


+150 


2007af 


+301 


2007sr 


+190 
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Fig. 5. — Number of SN la with a least TV spectra. The solid line 
is for all SN la, while the dotted line corresponds to SN la with no 
reliable estimate for the time of maximum. 
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Fig. 6. — Rcdshift distribution for the 462 SN la in our sample. 
The hatched histogram corresponds to the SN la sample published 
by Matheson et al. (2008). There are 4 SN la at rcdshifts larger 
than 0.075: SN 1999cc, SN 2006bu (part of the CfA3 light-curve 
sample), SN 2001bp, and SN 1996ab (part of the CfAl light-curve 
sample). 



SN 


Host 


z 


Rcf.^' 


1993ac 


CGCG 307-023 


0.0503'" 


abs 


1996ab 


Anonymous 


0.123'^ 


SNID 


1998ex 


MCG +11-10-16 


0.0365 


em 


2000dl 


UGC 1191 


0.0611 


abs 


2001eo 


UGC 3963 


0.066 


SNID 


2001CS 


2X J2020834+1905246 


0.0420 


em 


2002hu 


MCG +06-6-12 


0.0374 


em 


2003fa 


MCG +07-36-33 


0.0404'^ 


em 


2003hw 


2X J03014982+3544343 


0.0420 


abs 


2004gz 


MCG +10-23-45 


0.0137 


em 


2005ch 


Anonymous 


0.027 


SNID 


20051z 


UGC 1666 


0.040 


SNID 


2006bu 


2X J13524703+0518496 


0.081 


SNID 


2007H 


Anonymous 


0.0426 


cm 


2007cq 


2X J22144070+0504435 


0.0262 


cm 


2007if 


Anonymous 


0.0742° 


1 


2007kf 


Anonymous 


0.044 


SNID 


2007kg 


2MFGC 18005 


0.006 


SNID 


2007qe 


Anonymous 


0.0239 


cm 


2008E 


Anonymous 


0.034* 


2 


SNF20080522-000 


Anonymous 


0.0450S 


3 


SNF20080522-011 


Anonymous 


0.0379S 


3 


SNF20080623-001 


Anonymous 


0.0430S 


3 


SNF20080720-001 


Anonymous 


0.0209S 


3 



References. — (1) Scalzo ot al. 2010; (2) Yuan ot al. 2008; (3) 
Bailey ot al. 2009. 

Note. — 2X = 2MASX. 

^cni— emission lines from the host galaxy in our own speetruni; 
abs— absorption lines from the host galaxy in our own spectrum; 
SNID— cross-correlation with SN la spectral templates using the 
SNID code of Blondin & Tonry (2007) [typical error < 0.005]. 

NED reports a rcdshift z— 0.049 based on the host-galaxy redshift 
given in lAUC 5882 (Schmidt et al. 1993). We have re-derived the 
redshift using the same spectrum of the host galaxy as Schmidt et al. 
(1993). 

°NED reports an approximate redshift z — 0.13 leased on the 
blueshift of the Si II A5972 line (assumed to be lOOOOkms^"' in 
lAUC 6405; Garnavich et al. 1996); Riess et al. (1999) give log(cz) = 
4.571, or cz = 37239 kms"^ (z = 0.1242), "from [their] spectra of 
the host galaxies"; last, Jha et al. (2007) give cz — 37109 km s~ 
(z — 0.1238). There are no obvious galaxy lines in our spectrum of 
SN 1996ab, nor could we find a spectrum of its host galaxy. 

NED reports an erroneous redshift of cz — 1800km s^ (z — 
0.0060) for MCG +07-36-33. 

"Scalzo et al. (2010) report z = 0.07416 ± 0.00082 based on 
[O ll] A3727 and Ha. 

Rcdshift of a nearby galaxy cluster; the SN spectrum is consistent 
with this redshift. 

^Bailey et al. (2009) report the CMB-framo redshift for this SN. 
We derived a heliocentric redshift using its coordinates as reported 
on the SNfactory Supernova Discoveries webpage. 



1991T/1999aa-likc spectra, while 34 have spectra similar 
to SN 1991bg or SN 1986G. Of the three SN la with 
the smallest Ami5(i?) in our sample, two (SN 2006gz, 
Hicken et al. 2007b and SN 2007if, Scalzo et al. 2010) 
were speculated to originate from super-Chandrasekhar- 
mass progenitors. 

4. SPECTROSCOPIC SUBCLASSES 

We investigate how SN la can be classified into various 
subclasses based on their spectra, using the classification 
schemes of Branch et al. (2006) and Wang et al. (2009a). 
We also attempt to quantify the spectroscopic and pho- 
tometric variation within each subclass. Last, we present 
updated correlations of several spectroscopic indicators 
with A77ii5(i?). The results in this section are largely 
based on spectra from the CfA SN Program (see § 3), 
but we complement them with published data from the 



Diversity of SN la Spectra 



50 



40 



5 30 



i 20 



10 



Olay 
06gz 
07if 
051u 



\y^/i Matheson et al. (2008) 



0.5 




07hj 

1^ Ko6bz 

r7:7V//}/ i / / /., / /J / , 1 /7JA . I 



1.5 
Amj5(B) [mag] 



2.0 



Fig. 7. — Histogram of the 247 SN la for which we have a rcUable 
Amirj(B) measurement. The hatched histogram corresponds to 
the SN la sample published by Matheson et al. (2008). There 
arc four SN la with Ami5(B) ^ 0.7mag (SN 2001ay, SN 20051u, 
SN 2006gz, and SN 2007if) and three SN la with Ami5{B) ^ 2 
(SN 1999by, SN 2006bz, and SN 2007hj). 

literature, whose references are given in Table A2. 

4.1. Branch and Wang classification schemes 
4.1.1. Branch et al. (2006) classification 

Branch et al. (2006) presented a new classification 
scheme for SN la based on measurements of the pseudo 
equivalent widths (pEW) of absorption features near 
5750 A and 6100 A (generally attributed to Si iiA5972 
and SiiiA6355, respectively) in their maximum-light 
spectra. The overall shape of the Si ii A6355 absorption 
also enters this classification scheme. They divided SN la 
into four subgroups based on their position in this two- 
dimensional parameter space, but noted the absence of 
strict boundaries between them, suggesting that the dif- 
ferent subgroups were not physically distinct but instead 
corresponded to a continuous distribution of properties 
(with the possible exception of SN 2002cx-likc events) . 

Figure 8 {left) shows the pEW of the Si iiA5972 line 
vs. that of the Si ii A6355 line for 218 SN la (see Fig. 2 of 
Branch et al. 2009), based on spectra within 5 days from 
maximum light (see Table 4) . The different symbols cor- 
respond to the various spectroscopic subclasses defined 
by Branch et al. (2006). Example maximum-light spectra 
within each subclass are shown in Fig. 9. The exact loca- 
tions of the boundaries between the different subclasses 
are ill-defined, but this does not affect our analysis. 

The cluster of black dots forms the "Core Normal" 
(CN) subclass. SN la with similar pEW(Si ii A5972) 
but with broader Si ii A6355 absorptions (characteristic 
of larger expansion velocities) are labeled "Broad Line" 
(BL). Together, these subclasses constitute the bulk of 
SN la loosely referred to as "normal" in the literature. 
There is no strict boundary between the CN and BL 
subclasses, although the SN la we classify as BL have 
pEW(SiiiA6355)> 105 A. As noted by Branch et al. 
(2009), there is not a one- parameter sequence from the 
CN subclass to the more extreme BL SN la (such as 
SN 1984A, SN 2006X, and SN 2006bq). Two SN la clas- 
sified as BL by Branch et al. (2009) are re-classified by 
us as CN (SN 1999gd, SN 2002er; see Table 4) based on 
the magnitude of their pEW(Si ii A6355). We also clas- 
sify SN 2009ig in the CN subclass, as opposed to the BL 



subclass (see Parrent et al. 2011). 

The "Cool" (CL) subclass consists of SN la with 
deeper Si ii A5972 absorptions, often associated with low- 
luminosity events. They are referred to as "Cool" due to 
the notable absorption band around ^ 4200 A caused by 
lines of Ti ii, a signature of lower temperatures in the 
line- forming region (see, e.g., Hatano et al. 1999). The 
strength of this absorption band varies greatly within 
the CL subclass: it is strong in SN 2005ke, moderate 
in SN 1986G, and close to nonexistent in SN 2007au. 
The CL subclass contains SN 2006bt, which has a slowly- 
declining light curve characteristic of luminous SN la but 
with spectra displaying Ti ii absorption features charac- 
teristic of low- luminosity SN la (Foley et al. 2010). It 
also includes SN 2004eo, labeled "transitional" by Pa- 
storello et al. (2007) due to its intermediate properties 
between normal and low-luminosity SN la. Based on 
its resemblance to SN 2004eo, Branch et al. (2009) clas- 
sified SN 1989B as a CL SN la despite its shallower 
Si II A5972 absorption. This SN also displays narrower 
features in the 4600-5100 A spectral region, indicative of 
lower ejecta expansion velocities. However, the width 
of its Si II A6355 absorption is comparable to other BL 
SN la, and we choose to include SN 1989B in the BL sub- 
class (the developer of this classification scheme concurs; 
Branch 2011, priv. comm.). 

As seen from the inset in Fig. 8 (left:), there appears to 
be a continuous distribution between CL SN la at the low 
end of the pEW(Si ii A5972) distribution and CN and BL 
SN la at the high end of their pEW(Si ii A5972) distri- 
butions. However, we identify one CL SN la (SN 2007al) 
with a remarkably shallow silicon absorption. We include 
this SN in the CL subclass based on the presence of a 
prominent Ti ll absorption feature in its maximum-light 
spectrum (see Fig. 9). 

The last subclass, "Shallow Silicon" (SS), con- 
sists of SN la with weak SiiiA5972 and Si iiA6355. 
It is a heterogeneous category, including luminous 
1991T/1999aa-like events characterized by higher ion- 
ization lines (Fe ill), SN la resulting from possible 
super-Chandrasekhar-mass progenitors (SNLS-03D3bb 
or SN 2003fg, SN 2006gz, SN 2007if, SN 2009dc), faint 
2002cx-like SN (SN 2005hk, SN 2008A), and otherwise 
spectroscopically "normal" events (e.g., SN 2006S). Sev- 
eral SN la classified as SS by Branch et al. (2009) are re- 
classified by us as CN (SN 1999ac, SN 1999ec, SN 2000E, 
SN 2005cg; see Table 4) based on the magnitude of their 
pEW(Si IIA6355). SN 1999ac was ahcady noted as a 
borderline SS/CN by Branch et al. (2007), and its pecu- 
liar nature discussed at length by Garavini et al. (2005) 
and Phillips et al. (2006). 

Of the 246 SN la to which we assign a Branch et al. 
(2006) classification in Table 4, 94 (38.2%) are of the CN 
subclass, 74 (30.1%) are of the BL subclass, 43 (17.5%) 
are of the CL subclass, and 35 (14.2%) are of the SS 
subclass. Note that these numbers simply reflect the 
properties of the SN sample studied in this paper and 
not of the SN la class as a whole. 

By ^ 3 weeks past maximum light, spectra of SS SN la 
arc not distinguishable from those of the CN subclass at 
similar phases. By comparing the +19 d spectrum of the 
SN 1999aa to spectra of Core Normal SN la. Branch et al. 
(2009) found the best match to be SN 1994D at 4-14 d. 
The corresponding phase ratio (19/14 « 1.36) closely 
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Fig. 8. — Left: Pseudo-EW of the Si llA5972 absorption vs. that of the Si II A6355 absorption within five days from maximum Ught. 
When several measurements are available for a given SN within this phase range, we only show the measurement closest to maximum light. 
The different symbols correspond to the various spectroscopic subclasses defined by Branch et al. (2006). The inset corresponds to the 
region bounded by the dotted box in the main plot. Right: Pseudo-EW of the Si II A6355 absorption vs. the absorption velocity of the 
same line within three days from maximum light. When several measurements arc available for a given SN within this phase range, we 
only show the measurement closest to maximum light. The different symbols correspond to the various spectroscopic subclasses defined by 
Wang et al. (2009a). Note that peculiar SN la (e.g., SN 2002cx, SN 2005cc, SN 2005hk, SN 2008A) are not included in this classification 
scheme. The dotted lines correspond to the mean pseudo-EW and mean 11^1,3 of the Si II A6355 line for the Normal sample {open circles). 



TABLE 4 

SN lA CLASSIFICATION USING THE SCHEMES OF BRANCH ET AL. (2006) AND WANG 

ET AL. (2009a) 



SN 


pEW(A5972)=' 


pEW(A6355)'' 


i.abs(A6355)^ 


Phase'i 


Branch" 


Wangf 




(A) 




(A) 


(kms-i) 


(d) 


class 


class 


1981B 




9.1 


128.6 


-12050 


-1.5 


BL 


N 


1984A 




20.4 


200.5 


-16251 


-2.8 


BL 


HV 


1986G 




44.4 


122.7 


-10128 


-0.2 


CL 


91bg 


1989Bt 




20.9 


120.6 


-10634 


-1.3 


BL 


N 


1990N 




10.9 


87.1 


-9352 


1.6 


CN 


N 


1990O 




6.7 


91.1 


-12069 


-0.9 


CN 


N 


1991M 




18.4 


134.4 


-12458 


2.7 


BL 


HV 


1991T 




0.9 


29.0 


-9660 


-1.5 


SS 


91T 


1991bg 




47.9 


95.0 


-10004 


0.9 


CL 


91bg 


1992A 




19.1 


108.9 


-14192 


-0.8 


BL 


HV 



Note. — Tabic 4 is published in its entirety in the electronic edition of The Astronomical 
Journal. A portion is shown here for guidance regarding its form and content. 

^Pseudo-EW of the Si llA5972 absorption feature (formal uncertainty typically < lA). 
Pscudo-EW of the Si llA6355 absorption feature (formal uncertainty typically < lA). 

"^Absorption velocity the Si ll A6355 absorption feature (formal uncertainty typically < 
100 kms"^). 

Phase of measurements in rest-frame days from B-band maximum. 

■^Branch class based on pEW(A5972) and pEW(A6355) within 5d from maximum light: 
CN— Core Normal, BL— Broad Line, SS— Shallow Silicon, CL— Cool. We also report a classi- 
fication for SN la with no spectra within 5d from maximum (but within 7d) based on the 
overall shape of the Si II A6355 absorption feature. 

Wang class based on ^abs(A6355) within 7d from maximum light (note that Fig. 8 only dis- 
plays measurements within 3d from maximum): N— Normal, HV— High-velocity, 91T— 1991T- 
like, 91bg^l991bg-like. Peculiar SN la (SN 2000cx, SN 2002cx, SN 2003fg, SN 2005hk, 
SN 2006gz, SN 2007if, SN 2008A, SN 2008ae, SN 2009dc) are not part of this classification 
scheme. 

^SN 2003fg is also known as SNLS-03D3bb (Howell et al. 2006). 

'Reddening and contamination by the host galaxy renders the Branch classification uncer- 
tain. 

"''Branch classification differs from that reported by Branch et al. (2009) [indicated in paren- 
theses]: SN 1989B^BL (CL), SN 1999ac^CN (SS), SN 1999ee^CN (SS), SN 1999gd^CN 
(BL), SN 2000E^CN (SS), SN 2002er^CN (BL), SN 2005cg^CN (SS). 
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Fig. 9. — Montage of representative SN la speetra at niaxinium 
light in each of the spectroscopic subclasses defined by Branch 
et al. (2006). We indicate tentative line identifications correspond- 
ing to the dominant ion contribution to specific absorption fea- 
tures (the BL subclass shares the same line identifications as the 
CN subclass) . The vertical dotted line corresponds to Si II A6355 
blueshifted by 10000 km s~. All the spectra are from the CfA 
SN Program, except for the spectrum of SN 2009dc (Taubenberger 
et al. 2011). 



matched the ratio of their B-band hght curve stretch pa- 
rameters (1.143/0.838 w 1.36), prompting Branch et al. 
(2009) to suggest that SN la from the CN and SS sub- 
classes "age spectroscopically at the same rate that they 
decline photometrically." We have tested this hypothesis 
by finding the best-match CN SN la template to spec- 
tra of SS SN la in the phase range [-1-19, +23] d using 
SNID, and plot the ratio of their phases against that of 
their stretch parameters (from SALT2) in Fig. 10. While 
most SS SN la indeed lie on or close to the 1:1 relation, 



u 

C/5 






1.4 



1.2 



1.0 



0.8 - 



0.6 



08Z/03H 



'-V 



05eu/05bn^- 02hu/05na 
♦'06S/01fe 

,99gp/01cp 

*91T/07F 
'OlV/Olcp 



0.6 0.8 1.0 1.2 

Stretch ratio (SS/CN) 



1.4 



Fig. 10. — Ratio of phases of SN la spectra at three weeks 
past maximum in the SS subclass to that of the best-match SNID 
template in the CN subclass vs. the ratio of their _B-band light- 
curve stretch parameters from SALT2. The labels associated with 
each point indicate the name of the SS supernova followed by the 
name of the best-match CN template. The dashed line is the 1:1 
correspondence. 

there are a few significant outliers (including the peculiar 
SN 2000CX and the 1991T-like SN 1998ab), and the Pear- 
son correlation coefficient between both ratios is only 
r = 0.36. We find similar results when considering SS 
SN la at two and four weeks past maximum light, respec- 
tively. The hypothesis of Branch et al. (2009) therefore 
does not apply universally to SN la of the SS subclass. 

The spectroscopic properties of SN la belonging to 
these different subclasses have been discussed at length 
by Branch et al. (2006, 2007, 2008) for maximum-light, 
pre-maximum, and post-maximum epochs, respectively, 
and by Branch et al. (2009) at all epochs based on a 
larger sample. Here we focus on the spectroscopic vari- 
ation at maximum-light. Within each subclass, we gen- 
erate composite spectra using the same pre-processing 
as done by SNID: the individual spectra are "flattened" 
through division by a pseudo continuum. We then com- 
pute the mean fiux in each wavelength bin, as well as 
the standard (and maximum) deviation from the mean. 
The result is a composite spectrum with error bands for 
each subclass, which we show in Fig. 11. The large vari- 
ation within the SS subclass is clearly visible through- 
out the optical range (3500-7500 A). The CN subclass 
displays the smallest variation at any given wavelength, 
as expected, albeit with a noticeable variation in the 
strength of the Ca ii H&K absorption feature. The com- 
posite spectrum for the BL subclass reveals the presence 
of high-velocity components to the Si ii A6355 line in the 
maximum-deviation spectrum (dark blue) in addition to 
the large variation blueward of ^ 4000 A as for the CN 
subclass (see, e.g, SN 2001bf and SN 2007jg in Fig. 9). 
The CL composite spectrum shows the varying strength 
of the Ti II absorption feature as well as the large range in 
depth and position of the Si ii A6355 absorption (shallow 
in SN 2007al, less blueshifted in SN 2002es; see Fig. 9). 

In a recent paper, Stritzinger et al. (2011) noted a sim- 
ilarity between the photometric properties of SN 2006ot 
and SN 2006bt (Foley et al. 2010): a broad light curve 
characteristic of luminous SN la (Ami5(i?) = 0.84 mag), 
but a weak secondary i-band maximum characteristic 
of low-luminosity events. They also noted major differ- 
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Fig. 11. — Composite maximum-light spectra for the various 
spectroscopic subclasses defined by Branch et al. (2006). Individ- 
ual spectra have been flattened through division by a pseudo con- 
tinuum (see text for details). The shaded bands correspond to the 
standard {light blue) and maximum [dark blue) deviation about the 
mean maximum-light spectrum. The number in between parenthe- 
ses indicates the number of individual spectra that were used to 
generate each composite spectrum. 

ences between both SN in their maximum-hght spectra. 
SN 2006bt has a deep Si ii A5972 absorption and evidence 
for absorption by Ti II, both of which are characteristic 
of low-luminosity SN la. SN 2006ot, on the other hand, 
is characterized by a shallower Si ii A5972 absorption and 
a broad Si ii A6355 absorption, with no evidence for fea- 
tures associated with Ti ii. In the classification scheme 
of Branch et al. (2006), SN 2006bt belongs to the CL sub- 
class while SN 2006ot is an extreme BL SN la, reflecting 
these spectroscopic differences. Nonetheless, Stritzinger 
et al. (2011) argue that at 3-4 weeks past maximum light, 
the spectra of SN 2006ot are more similar to those of 
SN 2006bt than to those of the normal SN 2006ax (part 



of the CN subclass), again stressing the similarity be- 
tween both SN. At such late phases, however, differences 
between the various subclasses are less apparent than at 
maximum light, as noted by Foley et al. (2010)^^. We 
have run SNID on the +24 d spectrum of SN 2006ot and 
find good matches to spectra of SN la from the CN, CL, 
and BL subclasses at similar phases. The best-match 
spectrum of SN 2006bt (at -h33d) is ranked only 21"* 
overall. Similarly, the best-match template spectrum for 
the +33 d spectrum of the CL SN 2006bt (Foley et al. 
2010) is the CN SN 2004S (Krisciunas et al. 2007), illus- 
trating the difficulty in distinguishing between different 
SN la subclasses at late times. If anything, SN 2006bt 
and SN 2006ot demonstrate how two SN la with similar 
light-curve shapes can have very different maximum-light 
spectra. 

4.1.2. Wang et al. (2009a) classification 

The classification scheme of Wang et al. (2009a) is 
based on the velocity at maximum absorption (or absorp- 
tion velocity, Wabs) of the Si iiA6355 line around maxi- 
mum light. It was designed to study the spectroscopic 
diversity amongst "normal" SN la. We use our own mea- 
surements on the same sample of 10 SN la used by Wang 
et al. (2009a) to define a mean Wabs trend within one 
week from B-band maximum. SN la with Si ii A6355 
absorption velocities 3cr above this mean trend (in abso- 
lute terms) are classified as "High- velocity" (HV), the 
remaining objects forming the "Normal" (N) subclass 
(see Table 4). Note that this classification scheme does 
not apply to luminous 1991T/1999aa-like SN la, nor 
to low-luminosity 1991bg-like events or peculiar objects 
(2002cx-like, possibly super- Chandrasekhar etc.). Of the 
246 SN la to which we assign a Wang et al. (2009a) clas- 
sification in Table 4, 143 (58.1%) are of the N subclass, 
while 51 (20.7%) belong to the HV subclass. Again, these 
numbers simply reflect the properties of the SN sample 
studied in this paper and are not meant to reflect the 
SN la class as a whole. 

This classification scheme is illustrated in Fig. 8 (right), 
where we show the relation between the pseudo-EW of 
the Si II A6355 line vs. its absorption velocity for SN la 
with spectra within 3 d from maximum light. This figure 
is largely similar to Fig. 2 of Wang et al. (2009a), the Nor- 
mal SN la forming a cluster in the (pEW,t'abs) plane, and 
a tail of High- velocity objects following a trend of larger 
pEW for larger absorption blueshifts. The mean pEW 
and mean Uabs for the Normal sample are ~ 95 A and 
~ —11000 km s~ , respectively, comparable to the values 
found by Wang et al. (2009a). At maximum light, the 
boundary between the N and HV subclasses is located 
at ~ —12200 km s~ . There is a large overlap between 
the Core Normal subclass of Branch et al. (2006) and 
the Normal subclass of Wang et al. (2009a), as well as 
between the BL and HV subclasses (with a few excep- 
tions such as SN 2009ig, which is part of the CN and HV 
subclasses). The 1991T-like objects have absorption ve- 
locities similar to the Normal subclass, but with smaller 
pEW. AU 1991T-like SN la belong to the Shallow Sflicon 
subclass of Branch et al. (2006), except for the 1999aa- 

^^ In fact, Foley et al. (2010) note the similarity between the 
spectra of SN 2006bt and SN 2006ax at ~ 40 d past maximum 
light. 
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like SN 1999ac. SN la with spectra similar to SN 1991bg 
also have similar velocities but on average larger pEW. 
Notable exceptions include SN 2007al with similar ve- 
locities but much smaller pEW (also an outlier in the 
CL subclass; see Fig. 8, left), and SN 2002es with similar 
pEW but almost a factor of two lower |wabs|- This latter 
SN is discussed in more detail by Ganeshalingam et al. 
(2012). Ah 1991bg-like SN la belong to the Cool subclass 
of Branch et al. (2006). 

The largest measured Si ii A6355 pEW and blueshift 
in Fig. 8 is for SN 1984A (pEW« 200 A, |wabs| « 
16250 km s^ ). This SN is also the most extreme BL 
SN la in the classification scheme of Branch et al. (2006). 
We note that Wang et al. (2009a) measure pEW« 170 A 
and |wabs| ~ 15200kms~^ (see their Fig. 2), but their 
measurement is based on a non-flux calibrated spectrum 
directly scanned from Fig. 4 of Barbon et al. (1989), bi- 
asing their pEW measurement to lower values. More- 
over, their use of the classical Doppler formula (as op- 
posed the rclativistic version used in this paper; see Eq. 6 
in Blondin et al. 2006) leads to a ~ 500kms~^ lower 
blueshift when converting the wavelength at maximum 
absorption in Si ii A6355 to a velocity. 

4.2. Photometric properties as a function of 
spectroscopic class 

Based on a restricted sample of 9 SN la from the CN 
subclass. Branch et al. (2009) inferred a mean absolute 
peak i3-band magnitude of —19.48 mag with a scatter of 
only 0.14 mag, suggesting that SN la belonging to this 
subclass could be true standard candles. They noted, 
however, that this should be tested on a larger sample of 
SN la in the Hubble flow. 

We use the BayeSN statistical models of Mandcl et al. 
(2009, 2011) to infer the intrinsic peak absolute S-band 
magnitude (and hence the total extinction along the line 
of sight) of SN la in a consistent manner (assuming 
a Hubble constant of Hq = 72kms^^ Mpc^^; see Ta- 
ble 1). These statistical models describe the apparent 
distribution of light curves as a convolution of intrin- 
sic SN la variations and a dust distribution. Mandel 
et al. (2011) modeled the intrinsic covariance structure 
of the full multi-band light curves, capturing population 
correlations between the intrinsic absolute magnitudes, 
intrinsic colors and light curve decline rates over multi- 
ple phases and wavelengths, as well as the distribution 
of host galaxy dust and an apparent correlation between 
the dust extinction Ay and its wavelength dependence, 
parameterized by Ry- The models fit individual optical 
and NIR SN la light-curve data to estimate the dust ex- 
tinction, apparent and absolute light curves, and intrin- 
sic colors for each SN. These models were trained on a 
nearby (z < 0.07) set of SN la with optical (CfA3, Hicken 
et al. 2009; Carnegie SN Program, Contreras et al. 2010) 
and NIR (PAIRITEL; Wood-Vasey et al. 2008) data, plus 
light curves from the literature with joint optical and 
NIR observations. For the fits used in this paper, we 
have employed the distance-redshift constraint, except 
where alternative distance information was used (e.g., 
Cepheids; see § 7.2). 

Figure 12 shows the resulting intrinsic peak Mb dis- 
tribution for SN la at CMB-frame redshifts zcmb > 0.01 
with an inferred visual extinction Ay < 1 rnag, for the 
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Fig. 12. — Normalized distributions of intrinsic peak absolute 
B-band magnitude for the various spectroscopic subclasses defined 
by Branch et al. (2006) [top] and Wang et al. (2009a) [bottom]. The 
peak rest-frame magnitude and host-galaxy extinction were deter- 
mined using the BayeSN light-curve statistical model of Mandel 
et al. (2011). Only SN la with a CMB-frame redshift ^CMB > 0.01 
and an inferred visual extinction Ay < 1 mag were included. In 
each panel we indicate the sample mean and standard deviation of 
the intrinsic peak Mg (sec also Table 5). The number of SN la in 
each subclass is indicated in between parentheses. We do not show 
the distributions for the CL or 91bg subclasses since the BayeSN 
statistical model did not include SN la with Ami5(B) > 1.6 mag 
in its training set. The same is true of faint 2002cx-like SN la 
and luminous (possibly supcr-Chandrasekhar) events from the SS 
subclass. 

various spectroscopic subclasses defined by Branch et al. 
(2006) [top] and Wang et al. (2009a) [bottom]. The 41 
SN la in the CN subclass that satisfy these criteria have 
a mean intrinsic peak Mb = —19.40 mag with a{MB) = 
0.16 mag. This is comparable to the values derived for 
the Normal subclass of Wang et al. (2009a) [Mb = 
— 19.38 mag; u^Mb) = 0.17 mag]. The intrinsic peak Mb 
scatter is small for the CN (or Normal) subclass, but it 
is in fact larger than for the SS {a{MB) = 0.13mag) 
and BL subclasses (i7(Ms) = 0.15 mag) [the same is also 
true when comparing the Normal subclass to the 91T and 
HV subclasses of Wang et al. 2009a]. We note, however, 
that the SS distribution does not include faint 2002cx- 
like SN la (which typically have Mb ^ —18 mag; e.g., 
Mb « -17.5 mag for SN 2002cx and Mb « -18.0mag 
for SN 2005hk; see Phillips et al. 2007) nor luminous 
(possibly super-Chandrasekhar) events (e.g., SN 2006gz 
with Mb « -19.9, Hicken et al. 2007a; or SN 2009dc 
with Mb ~ —20.2, Taubenberger et al. 2011) since 
BayeSN did not include such SN la in its training set. 

The mean intrinsic peak Mb increases steadily along 
the SS^-CN— >-BL sequence, with a corresponding in- 
crease in the mean A?Tii5(i?) [see Table 5]. The same 
applies to the 91T— >N— i-HV sequence. There is a hint of 
redder mean intrinsic B — V color for the BL (HV) sub- 
class compared to the CN (N) subclass, as expected from 
the "brighter-bluer" relation of Tripp (1998), but the dif- 
ference (0.02-0.03 mag) is small compared to the scatter 
(0.06-0.07mag) within each subclass and only measur- 
able with sample averages. Low-luminosity SN la from 
the CL or 91bg subclasses are not included in this com- 
parison since the BayeSN statistical model did not in- 
clude SN la with Ami5{B) > 1.6 mag in its training set. 
Supcrnovae similar to SN 1991bg have intrinsic peak B- 
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TABLE 5 

Mean properties of SN Ia in the classification schemes of 

Branch et al. (2006) and Wang et al. (2009a) 



Class 


(A/b> ((t) 


(Ami5(iJ)> (<7) 


{B - V) (a) 


A^SN 


CN 
BL 


-19.57 (0.14) 
-19.40 (0.16) 
-19.30 (0.15) 


0.94 (0.13) 
1.12 (0.18) 
1.22 (0.15) 


-0.11 (0.07) 
-0.12 (0.06) 
-0.10 (0.07) 


21 
41 
37 


91T 

Normal 

HV 


-19.59 (0.13) 
-19.38 (0.17) 
-19.28 (0.16) 


0.93 (0.14) 
1.15 (0.21) 
1.21 (0.14) 


-0.12 (0.07) 
-0.12 (0.06) 
-0.09 (0.07) 


13 
67 
23 



Note. — SS=Shallow Silicon; CN=Coro Normal; BL=Broad Line; 
HV— High-velocity. Cool (CL) and 1991bg-like SN la are not shown 
here (see text for details). 

^excluding 2002cx-like SN la and super-Chandrasekhar events. 



band magnitudes ~ 2 mag fainter than normal SN la, 
with A77ii5(i?) > 1.9 mag, and much redder intrinsic 
B — V colors at maximum light. A prime example is 
SN 2005bl, with Mb « -17.2 mag, Ami5(S) k. 1.9 mag, 
and [B — V^)max ~ 0.5 mag (Taubenberger et al. 2008). 

One might wonder whether SN la in different sub- 
classes obey different width-luminosity relations (WLR) . 
Figure 13 [left] shows the intrinsic peak Mb inferred 
from BayeSN fits vs. Amir,{B) for the SS, CN, and BL 
subclasses of Branch et al. (2006). Overplotted are best- 
fit relations of the form Mb — a [Amis (B) — 1.1] + b for 
the entire sample (solid line) and for the individual sub- 
classes [dashed and dotted lines). The coefficients of the 
linear fits are given in each case. The slope of the WLR 
gets steeper along the SS— ?>CN— >BL sequence, while the 
intercept corresponds to a progressively fainter fiducial 
intrinsic peak Mb along this same sequence. The slopes 
and intercepts for the individual subclasses are however 
consistent within ^ la of one another, such that a single 
average WLR is adequate to describe all the subclasses 
(other than the CL subclass). As seen from the right 
panel of Fig. 13, the same analysis holds when replac- 
ing (SS,CN,BL) with (91T,N,HV) for the Wang et al. 
(2009a) classification scheme, the only difference being 
the intercept for the HV subclass which is ~ 2cr larger 
than those for the Normal and 91T subclasses. 

4.3. Spectroscopic luminosity indicators 

Several authors have studied the correlation between 
various spectroscopic indicators and light-curve param- 
eters, with the aim to use such indicators to improve 
distance measurements to SN la. Blondin et al. (2011b) 
showed that spectroscopic indicators alone could com- 
pete with the standard light-curve parameters for dis- 
tance measurements, but that combining spectra with 
photometry yielded no statistically significant improve- 
ment (see also Silverman et al. 2012b). Their analysis 
was based on a small number of supernovae, however, so 
that larger samples may reveal spectroscopic measure- 
ments that provide independent information on the lu- 
minosity of SN la. Such correlations are also potentially 
useful to constrain SN la models (see, e.g., Blondin et al. 
2011a), and we investigate the relation of three spectro- 
scopic indicators with Ami5(i?) in Fig. 14. 

The upper-left panel shows the 7?.(Si) ratio of Nugent 
et al. (1995) (defined as the ratio of the relative absorp- 
tion depth of the Si ii A5972 line to that of Si ii A6355), 



which is thought to define a temperature sequence from 
low-luminosity 1991bg-like SN la to luminous 1991T-like 
events (see also Hachinger et al. 2008). The different 
symbols correspond to the spectroscopic subclasses de- 
fined by Branch et al. (2006), as in the left panel of Fig. 8. 
There is a clear correlation of 7?,(Si) with Ami5(_B) (the 
Pearson correlation coefficient is r = 0.82), with some 
notable outliers including two SN la from the CL sub- 
class: SN 2002fb and SN 2006bt (whose peculiar nature 
has been discussed at length by Foley et al. 2010). At low 
Ami5{B), the BL SN 2001ay (the most slowly-declining 
SN la to date; Krisciunas et al. 2011b) and the pecu- 
har SS SN 2006gz (possibly resulting from a double WD 
merger; Hicken et al. 2007a) are modest outliers. 

The 7?.(Si) ratio has a relatively large associated mea- 
surement error, since the relative absorption depth is 
measured at a specific wavelength. Using the ratio of 
pseudo-EW of both lines instead not only significantly re- 
duces the error on individual measurements (since we are 
integrating over many wavelength bins), it also leads to 
a stronger and tighter relation with Ami^{B) {r — 0.86; 
upper-right). SN 2006bt remains an outlier using the 
pseudo-EW ratio, but SN 2002fb now follows the linear 
relation defined by the bulk of the sample. As noted em- 
pirically by Hachinger et al. (2006), and later confirmed 
theoretically by Hachinger et al. 2008, it is the Si ii A5972 
line that drives the correlation of 7^(Si) and its pseudo- 
EW equivalent with Ami5{B). Using the relative ab- 
sorption depth or the pseudo-EW of the Si ii A5972 line 
alone leads to equally strong correlations with Ami^{B) 
(r = 0.79 and r = 0.85, respectively; not shown). 

The lower-left panel shows the pseudo-EW of the 
Si IIA4130 line vs. Ami5(B) (see Bronder et al. 2008; 
Arsenijevic et al. 2008; Walker et al. 2011; Blondin et al. 
2011b; Chotard et al. 2011). The correlation is very clear, 
although a non-linear relation may need to be invoked 
to accommodate the CL SN la at Ami^{B) > 1.5 mag. 
The peculiar SN 2006gz is again a modest outlier, but 
SN 2001ay is now a significant outlier in this relation. 
The largest outlier however is the BL SN 2006ot (which 
was not a significant outlier in the TZ{Si) vs. Ami5{B) 
relation) , whose similarity in terms of photometric prop- 
erties with SN 2006bt was noted by Stritzinger et al. 
(2011). We argued in § 4.1.1 that both objects are spec- 
troscopically distinct, and this is clearly visible in this 
plot, where SN 2006bt follows the same relation as the 
bulk of the SN la sample. 

Finally, the lower-right panel shows the 7?.(Si, Fe) ra- 
tio of Hachinger et al. (2006), defined as the ratio 
of the pseudo-EW of the Si ii A5972 line to that of 
the Fe ilA4800 feature (see Blondin et al. 2011b, their 
Fig. 14). The correlation with Ami5{B) is again very 
strong (r = 0.80). In contrast to pEW(Si ii A4130), 
SN 2006bt is this time the largest outlier, while 
SN 2006ot lies close to the mean relation. As was the 
case for the 7^(Si) ratio, both SN 2001ay and SN 2006gz 
are modest outliers. Since 7?,(Si, Fe) and the pseudo-EW 
equivalent of 7^(Si) are both positively correlated with 
Ami5(i?), one expects the pseudo-EWs of FeiiA4800 
and Si ii A6355 (at maximum light) to also be corre- 
lated. There is indeed a modest correlation (r = 0.60) 
between pEW(Si ii A6355) and pEW(Fe ii A4800), but 
the scatter is large and the relation fiattens off for 
pEW(Fe iiA4800)> 150 A. 
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5. LINE VELOCITIES 

We now focus on the defining Si ii A6355 line of Type 
la supernovae, studying the variation in the time evo- 
lution of its absorption velocity (§ 5.1). We propose 
an alternative (and in our view more robust) method 
to Benetti et al. (2005) for measuring velocity gradients 
(§ 5.2). Last, we investigate the correlation of Si ii ve- 
locity with intrinsic color at maximum light suggested by 
Foley & Kasen (2011) in § 5.3. 

5.1. Evolution of the Si uX6355 absorption velocity 

Figure 15 shows the time evolution of the Si ii A6355 
absorption velocity for the different spectroscopic sub- 
classes of Branch et al. (2006). Blondin et al. (2006) 
presented similar measurements (in different ATOi5(i?) 
bins) along with a physical interpretation in terms of mi- 
croscopic atomic-transition properties and macroscopic 
properties of the ejecta (density and velocity distribu- 
tions). We refer the reader to that paper for more detail 
on the impact of these properties on the morphology of 
line profiles in SN la ejecta. 

For the Core Normal subclass {top left), there is a large 
variation in Dabs(Si ii A6355) at phases < — lOd (e.g., 
SN 2003du and SN 2006gr have v^^s « -12000 km s"^ 
and — 20000kms~^ at — lOd, respectively). By max- 
imum light, the spread in Wabs has decreased to < 
3000 km s~^, SN 2009ig being a clear outlier (it is part 
of the High-velocity subclass of Wang et al. (2009a), but 
the pseudo-EW of its Si ii A6355 line at maximum light 
places it unambiguously in the CN subclass of Branch 
et al. (2006); sec Fig. 8). The evolution is fast at early 
times and more gradual at post-maximum epochs, re- 
flecting the recession of the line-forming region in the co- 
moving frame of the ejecta. There is also a large spread 
in Vahs at early times for SN la in the Broad Line sub- 
class {top right), but unlike the CN subclass it remains 
large at later phases (> 5000 km s~^ variation at maxi- 
mum light). At phases later than — lOd, the absorption 
velocity is on average larger at any given time when com- 
pared to the CN subclass, illustrating the large overlap 
between the BL subclass and the High- velocity subclass 
defined by Wang et al. (2009a). 

The bulk of the Shallow Silicon subclass {bottom left) 
has velocities comparable to the CN subclass, except at 
early times where there are no measurements |uabs| > 
15000 km s"\ A handful of SN la display -5000 > 
I'abs j^ —10000 km s~ at maximum light. These include 
2002cx-like SN la (SN 2005hk, Wabs.max ~ -5500 km s"^; 
SN 2008A, Wabs,max « -7500kms"i), SN la re- 
sulting from possible super-Chandrasekhar progenitors 
(SN 2009dc, Uabs,max ~ -8000 km s'^; SN 2007if, 
Wabs.max « -7500 km s"^), and the 1991T-like SN 2005M 
(wabs.max ~ —8000 km s^ ). This last SN stands out 
amongst other 1991T-like SN la which all display > 
10000 km s^ absorption blueshifts at maximum light. 

The Cool subclass {bottom right) displays velocities 
similar to the CN and SS subclasses around maximum 
light, but with a more rapid evolution at post- maximum 
epochs. One significant outlier is SN 2002es (Gane- 
shalingam et al. 2012), with Wabs ~ — 6000kms~^ at +2 d 
(see also Fig. 8, right), comparable to the SS SN 2005hk 



at a similar phase. 

Benetti et al. (2005) noted the large variation in post- 
maximum evolution of the Si ii A6355 absorption velocity 
for different SN la. This is clearly visible in Fig. 15, not 
only across the different subclasses but also within each 
subclass (the largest variation being apparent within the 
BL subclass). Benetti et al. (2005) introduced the ve- 
locity gradient, noted v, representing the average \vi,hs\ 
decline rate past maximum light (units of km s~ day" ) , 
and divided their SN la sample into three groups (FAINT 
for 1991bg-like SN la, HVG and LVG for SN la with high 
and low velocity gradients, respectively). They suggested 
that the LVG and HVG groups could consist of SN la 
with different properties of their outer ejecta, resulting 
from different degrees of mixing or from interaction with 
circumstellar material. A recent study by Maeda et al. 
(2010a) associates the observed diversity in the velocity 
gradients with viewing-angle effects in off-center delayed- 
detonation models. Blondin ct al. (2011a) showed that 
varying the criterion for deflagration-to-detonation tran- 
sition in explosions resulting from an isotropic distribu- 
tion of ignition points (i.e., not off-center explosions) re- 
sulted in a range of velocity gradients comparable to that 
observed by Benetti et al. (2005). 

The notion of velocity gradients has received much at- 
tention since the paper by Benetti et al. (2005), and there 
are different physical interpretations as to the origin of 
their variation. This however can only be properly ad- 
dressed if the actual parameter is deflned in a robust way, 
which we set out to do in the following section. 

5.2. Redefining velocity gradients 

According to Benetti et al. (2005), the velocity gradient 
is derived "from a least-squares fit of the measurements 
taken between maximum and either the time the Si ii fea- 
ture disappears or the last available spectrum, whichever 
comes first." This definition makes two implicit (and 
related) assumptions: (1) the evolution of iiabs is a lin- 
ear function of time at post-maximum epochs, and (2) 
the phase of the last T;abs measurement (or the range of 
available measurements) has little impact on the derived 
velocity gradient. 

Figure 16 shows that both assumptions can have a sig- 
nificant impact on the derived gradient. In the left panel 
we show the time evolution of Wabs (Si ii A6355) for the BL 
SN 2002bo {filled circles) and the CN SN 2005cf {open 
squares) . While the Wabs evolution for SN 2005cf appears 
to be almost linear at post-maximum epochs, that of 
SN 2002bo is clearly not, with a more rapid change in 
the interval [-1-10,-1-15] d. Moreover, at phases > -f-25d 
the Si II A6355 absorption profile is increasingly contam- 
inated by the emergence of lines from different ions {see 
inset), such that the iiabs measurements at these times 
are not reliable. The right panel of Fig. 16 shows the 
variation in the measured velocity gradient in the in- 
terval [-l-0,iiast] d as a function of the phase of the last 
available measurement, iiast- The velocity gradient for 
SN 2002bo is ^ 115 km s^ day" when measured over 
the interval [-1-0,-1-10] d, and ^ 165kms~^ day~^ when 
measured over the interval [-1-0, +20] d. Including the bi- 
ased Wabs measurements at > +25 d yields a velocity gra- 
dient ^ 120kms~"'^ day^^, comparable to that obtained 
for [+0,+10] d, but with a x^ for the linear fit that is 
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Fig. 15. — Evolution of the Si II A6355 absorption velocity with time for the four spectroscopic subclasses defined by Branch et al. (2006). 
The color-coding denotes individual SN la. For the CN and BL subclasses we have only included SN la for which we have at least 10 Vabs 
measurements prior to -|-30 d from maximum light, while for the SS and CL subclasses only SN la with at least 5 measurements are shown. 

sons why no spectra would liave been obtained at later 
times (poor weather, telescope problems, conflict with 
other observations etc.), none of which should affect the 
inferred intrinsic properties of SN la. 

We also note that the velocity gradient for SN 2006bt 
reported by Foley et al. (2010) is in error. We mea- 



two orders of magnitude larger. Benetti et al. (2005) 
report a velocity gradient of 110 ± 7kms~ day^ for 
this SN, due to the inclusion of measurements out to 
+30 d (Benetti 2008, priv. comm.). For SN 2005cf, the 
quasi-linear variation of Wabs with time results in a mi- 
nor variation of the velocity gradient with tiast , from 



65kins day 



^ for iiast ~ +10 d to ~ 50kms May ^ 



for iiast ~ +30 d. Wang et al. (2009b) report a velocity 
gradient of 38 ± 5kms~^ day~^ for this SN using mea- 
surements out to +30 d. 

The classification of SN 2002bo and SN 2005cf respec- 
tively in the HVG and LVG groups still holds, regard- 
less of iiast and whether or not the biased Wabs measure- 
ments of SN 2002bo are included. This is not always 
the case. For SN 2002er, Benetti et al. (2005) report a 
velocity gradient of 92 ± 5 km s~ day" , including Uabs 
measurements past -1-20 d when the Si iiA6355 absorp- 
tion profile is clearly contaminated by other lines. Using 
measurements in the phase interval [-1-0, +16] d we find 
81±3kms~ day" , consistent with their value, and con- 
firming the classification of SN 2002er as a HVG SN la. 
However, restricting the phase interval to [+0,+10] d, we 
infer a velocity gradient of 47±7kms~^ day~^, consistent 
with the LVG group. One can think of a variety of rea- 



sure a velocity gradient of 93 ± 16kms day over 
the time interval [+0,+16]d, while Foley et al. (2010) 
report 51kms~^day~^ over the same interval. Foley 
(2011, priv. comm.) has revised this measurement to 
96kms~ day" . According to Fig. 7 of Foley et al. 
(2010), this new measurement places SN 2006bt well into 
the HVG group, where Foley et al. (2010) originally in- 
cluded this SN in the LVG group (albeit close to the 
border with HVG SN la). 

A more robust definition of the velocity gradient should 
consider Wabs measurements over a fixed phase range, as 
done for instance by Foley et al. (2011) over the interval 
—6 < t < -l-lOd using linear functions to describe the 
Va.hs{t) curves. This still does not account for non-linear 
variations in absorption velocity. As an alternative, we 
suggest the measurement of the mean velocity decline 
rate Awabs/At over some fixed time interval [io,^i]: 
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In practice there are not always measurements at pre- 
cisely to and ti, so one may choose to either select the 
^abs measurement closest in time to to and fi (within 
a small time interval, e.g., ±2d), or to interpolate the 
iiabs measurements at the desired phases, again mak- 
ing sure the Wabs(t) data is well sampled around to 
and ti. Figure 17 (left) shows the mean velocity de- 
cline rate for the Si iiA6355 line for the time interval 
[+0,+10]d vs. AniiQ^B), where we choose the Uabs 
measurement closest in time and within 2 d from these 
time boundaries. The clustering into the three groups 
(FAINT, LVG, HVG) is less obvious than in Fig. 3 
of Benetti et al. (2005), but the points occupy similar 
regions of parameter space, with a smaller scatter in 
Aiiabs/Aij+o.+io] for Ami5(_B) > 1.5 mag. The Broad 
Line subclass displays the largest variation in mean ve- 
locity decline rate (confirming what was visually appar- 
ent in Fig. 15), ranging from ^ — 20kms^ day^ for 
SN 2006cj to +215kms"May"^ for SN 2006X. Two 
other SN have Awabs/At[+o,+io] < 0, SN 1991T and 
SN 2005M (both part of the SS subclass). A combination 
of optical-depth effects (e.g., a Si iii— >-Si ii recombination 
wave) and measurement error (the Si ii absorption is rel- 
atively weak and occurs on a steep pseudo-continuum) 
could explain why the absorption blueshift apparently 
increases over this phase range. Several objects have 
Ai)abs/At[+o +10] consistent with zero, including the CN 
SN 2006gr. This was also the case for SN 2005hj, for 
which Quimby et al. (2007) invoke a shell-like density 
structure (formed within the context of a pulsating de- 
layed detonation or WD- WD merger scenario) to explain 
the ^ lOd-long plateau in absorption velocity of the 
SiiiA6355 1ine. 

The CN and CL subclasses appear to form a sequence 
of linearly increasing Awabs/Afr+o,+io] with Ami^^B). 
Both quantities are indeed strongly correlated (r = 0.80) 
for these two subclasses. Within the BL subclass, there is 



a strong correlation (r ~ —0.80) between the Si iiA6355 
absorption velocity at maximum light and its subsequent 
mean decline over the interval [-1-0,-1-10] d. In other 
words, SN la with the largest absorption blueshifts at 
maximum light also display steeper velocity gradients 
(see also Benetti et al. 2005; Wang et al. 2009a; Foley 
et al. 2011). 

In the limit {ti — to) — i- 0, Eq. 1 yields an instanta- 
neous velocity decline rate at a given time, dt'abs/'^t|j=t(,. 
Wc determine this using second or third-order polyno- 
mial fits to the Va,hs{t) data, as shown in Fig. 16 {left; 
solid curves). The instantaneous velocity decline rates at 
t^+Od for Si II A6355 in SN 2002bo and SN 2005cf are 
shown as the error bars in the right panel of Fig. 16. Fig- 
ure 17 (right) shows the same quantity for a larger sample 
vs. Amiz{B). The different subclasses occupy similar re- 
gions of parameter space as when considering the mean 
velocity decline rate (not surprisingly, since the mean and 
instantaneous decline rates are strongly correlated), but 
reveals an even larger variation within the BL subclass, 
with instantaneous velocity decline rates at maximum 
light as large as ^ 300 km s~ day" for SN 2001ay and 
- 400kms"May-^ for SN 2003W. SN 2005hj has the 
lowest instantaneous velocity decline rate at maximum 
light (14 ± 5kms~^ day~^), in line with the results of 
Quimby et al. (2007). 

5.3. Intrinsic color and Si II velocity 

Based on the results by Wang et al. (2009a), Foley & 
Kascn (2011) noted that SN la in the High- velocity sub- 
class were intrinsically redder (in B"^^^ — V"^^^ pseudo- 
color^*^) than SN la in the Normal subclass. They pro- 
vided a simple physical explanation using radiative trans- 
fer calculations by Kasen & Plewa (2007) based on the 
detonating failed deflagration (DFD) model of Plewa 
(2007), ahhough Blondin et al. (2011a) showed that this 
relation did not hold in the 2D delayed-detonation mod- 
els of Kasen et al. (2009). Using the Si ii A6355 absorp- 

'^^ difference between the B magnitude at B-band maximum and 
the V magnitude at V-band maximum. 
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to ~ 9(7 in their paper. The scatter about the best-fit 
hnear relation [dashed line) is the same to within ^ 10% 
for both samples, in conflict with the resuhs of Foley 
ct al. (2011), who infer a ~ 50% larger scatter for the 
HV sample. 

The robustness of the relation between Si ii veloc- 
ity and intrinsic color remains questionable given the 
large uncertainty associated with the determination of 
the latter quantity. Foley et al. (2011) fit for a relation 
between light-curve shape corrected peak absolute My 
(not corrected for host-galaxy reddening) and observed 
]jmax_ymax pgeudo-color (their Fig. 15), based on Eq. 1 
of Wang et al. (2009a). They assume that the horizon- 
tal scatter about this relation is entirely due to intrinsic 
variations in B"^'^^ — \/™3,x pseudo-color, where it in fact 
also includes some intrinsic scatter about the luminosity- 
A77T,i5(i?) relation. It is thus unlikely that their "intrinsic 
color estimates" are simply offset from the true intrin- 
sic color, as postulated in their paper. Moreover, Foley 
et al. (2011) do not propagate the errors associated with 
their nuisance parameters (a, Mzp, Ry, and C in their 
Eq. 15) in the intrinsic color error bars, which artificially 
enhances the statistical significance of their correlation. 
The same limitations apply to the recent study of this 
relation for high-rcdshift SN la by Foley (2012), with the 
additional caveat that their analysis makes use of an indi- 
rect measure of the observed B"^^^ — ymax pseudo-color. 

BayeSN uses a more elaborate statistical treatment of 
the population correlations between the intrinsic abso- 
lute magnitudes (BVRI, including JH when available), 
intrinsic colors and light-curve shape, as well as the dis- 
tribution and properties of host galaxy dust. BayeSN 
therefore provides an estimate of the true intrinsic color, 
with realistic error bars reflecting the uncertainty in the 
various population correlations. It is then not surprising 
that our estimates for the intrinsic B"^'^^ — V'^'^^ pseudo- 
color are only weakly correlated (r = 0.21) with those of 
Foley et al. (2011). Likewise, the mean intrinsic colors we 
derive for both Normal and High- velocity samples are sig- 
nificantly bluer (by ^ 0.08 mag and ^ 0.12 mag, respec- 
tively, when using the intrinsic B"^^^ — y max pseudo-color 

significance when using the Baycsian approach to linear regression 
of Kelly (2007). 



tion velocity measurements presented in this paper, Fo- 
ley et al. (2011) quantified the relation between Si ii ve- 
locity and intrinsic B"^'^^ — V"^^'^ pseudo-color at max- 
imum light for Normal and High-velocity SN la with 
1.0 < Atoi5(_B) < 1.5 mag. They find that the slope of 
the relation between both quantities is non-zero at the 
~ 9(7 level using a linear weighted least-squares fit (but 
only at the ~ 3.5(7 level when using a Baycsian approach 
to linear regression; see Kelly 2007), with a Pearson cor- 
relation coefficient r = —0.39. 

We show the relation between the intrinsic B — V color 
at B-band maximum light, as derived from BayeSN fits 
to SN la light curves (see Table 1), and the absorption 
velocity of the Si iiA6355 line in Fig. 18 (left) for the 
Normal and High-velocity subclasses. We restrict the 
Ar7ii5(i?) range to 1.0 < Ami5(i?) < 1.5 mag, as done 
by Foley et al. (2011). Moreover we use their interpolated 
velocity measurements at maximum light (their Table 1). 
The right panel shows histograms of intrinsic B~V color 
at maximum light for the Normal (A^sn = 51; open) and 
High-velocity (A^sn = 22; filled) subsamplcs. Of the 65 
SN la used by Foley et al. (2011) in their analysis, 42 (27 
Normal and 15 High- velocity) are part of our sample. 

A linear weighted least-squares fit to the data in Fig. 18 
yields the following relation between intrinsic B — V color 
and Si ii velocity at maximum light; 

B-V = -0.26(0.05)-0.013(0.005)x[ugi/103kms^i] mag, 

(2) 
with a Pearson correlation coefficient r = —0.35 (see Ta- 
ble 6), consistent with that found by Foley et al. (2011). 
Taking into account the uncertainty in the intrinsic color 
estimates via a Monte Carlo simulation results in almost 
no correlation (r w —0.20). The slope of the best-fit 
linear relation is shallower than that found in their pa- 
per, and with a larger associated error that refiects the 
larger uncertainty we derive on intrinsic color estimates 
(typically 0.05 mag using BayeSN cf. 0.03 mag for the 
estimates reported by Foley et al. 2011). The statistical 
significance of the linear relation we find is also signifi- 
cantly lower than reported by Foley et al. (2011). Our 
fit rejects a non-zero slope at the ~ 3a level^^, compared 

^^ Unlike Foley et al. (2011), we obtain the same statistical 
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Fig. 18. — Left: Intrinsic B ~ V color at B-band maximum light, as derived from BayeSN fits to SN la light curves, vs. the absorption 
velocity of the Si II A6355 line, for SN la in the Normal and High-velocity subclasses with 1.0 < Amis{B) < 1.5mag. The Pearson 
correlation coefficient is r = —0.35 (r fa —0.20 when taking into account the uncertainties in intrinsic color). The dashed line is the result 
of a least-squares linear fit, with a slope of 0.013 ± 0.005 (~ 2a different from zero). We also indicate the weighted mean intrinsic color 
{{B — V)) for both samples. Right: Distribution of intrinsic B — V color at maximum light for the Normal and High- velocity subsamples. 
A K-S test shows that the two samples are discrepant at the p = 0.023 < 0.05 significance level. 



TABLE 6 
Intrinsic color vs. Si iiA6355 absorption velocity. 



Intrinsic color 


Slope 


r (color) jv (color) ^--ia 


A(color)Hv,iV 


K-S Prob. 


B-V 

omax T/max 


-0.013±0.005 
-0.011±0.005 


This paper (A^n = 51, Ahv = 22) 
(2.9(t) -0.35 -0.115±0.006 -0.085±0.011 
(2.3(t) -0.26 -0.101±0.007 -0.070±0.012 


0.030±0.013 (2.3cr) 
0.030±0.013 (2.3(t) 


0.023 
0.140 


omax yma:>i 


-0.033±0.004 


Foley et al. 2011 (An = 46, A^hv = 19)=' 
(S.Scr) -0.39 -0.022±0.005 0.047±0.008 


0.069±0.010 (7.0(t) 


0.003 



^Wc have used our own derived boundary between the N and HV subclasses at ~ — 12200 kms (sec § 4.1.2). 



from BayeSN) than those using the color estimates of Fo- 
ley ct al. (2011). The larger uncertainties derived from 
BayeSN fits reflect the difficulty in obtaining accurate 
estimates of intrinsic color, properly accounting for the 
population distributions (or variance) of intrinsic color 
and dust extinction. 

The difference in the weighted mean intrinsic B ~ V 
color between the High-velocity and Normal samples is 
0.030±0.013mag, corresponding to a ~ 2(7 difference (cf. 
~ 7(7 in Foley ct al. 2011). A Kolmogorov-Smirnov (K- 
S) test shows that the two samples are discrepant at the 
p = 0.023 significance level, which is less than the canon- 
ical cutoff of 0.05 for statistical significance^^, with the 
caveat that this statistic does not consider the (large) 
uncertainty in intrinsic color. This is still significantly 
larger (i.e., less significant) than the K-S probability in- 
ferred for the Foley et al. (2011) estimates of intrinsic 

^^ Note that a smaller p-value indicates a higher level of statis- 
tical significance for the K-S test. 



^max _ ^max pgeudo-color {p = 0.003). 

Using the intrinsic B — V color at B-band maximum as 
opposed to the iJ™'^'^ — ymax pseudo-color does not affect 
our results nor their interpretation, as both quantities are 
strongly correlated (r = 0.72; sec Tabic 1). The weighted 
mean difference between intrinsic B™^^ — V^™^^ pseudo- 
color and intrinsic B — V color at i?-band maximum is 
only 0.02 mag with an rms of 0.04 mag. Pseudo-colors are 
difficult to interpret since, for an individual supernova, 
the l/-band may reach its maximum brightness between 
^ Od and ~ 4d later than the i3-band (see, e.g., Blondin 
et al. 2011a, their Fig. 9). We prefer to use the true color 
at a given time. The resulting correlation between Si Ii 
absorption velocity and intrinsic B™^^ — V™'^^ pseudo- 
color is slightly weaker than when using the B — V color 
at i?-band maximum light [r — —0.26). While the High- 
velocity sample is found to have a redder average intrinsic 
^max_^max pgeudo-color than the Normal sample at the 
^ 2(7 level (i.e. the same significance than when using 
the B — V color), a K-S test shows that the two samples 
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are consistent at the p = 0.14 > 0.05 significance level. 

A larger sample of SN la with combined optical and 
NIR photometry should enable a more precise determi- 
nation of the intrinsic colors of SN la (see Mandel et al. 
2011) and a more reliable study of the relation with Si ii 
velocity. 

6. EARLY-TIME SPECTRA 

Spectra taken shortly (within ^ 10 d) past explosion 
probe the outer layers of the SN ejecta and hence the 
completeness of C/0 burning in the progenitor WD. The 
velocity interval over which intermediate-mass elements 
absorb (in particular Si and Ca, see § 6.1 and § 6.2) and 
the detection of lines from (unburnt) carbon (see § 6.3) 
both provide constraints on possible explosion mecha- 



6.1. Maximum velocity of Sill absorption 

In § 5.1 wc investigated the time evolution of the 
Si II A6355 absorption velocity as a function of spectro- 
scopic subclass (see Fig. 15). This velocity corresponds 
to the location of maxim,um, absorption. Here we are in- 
terested in the maximum velocity at which Si ii absorbs, 
since this value gives an estimate of the extent of incom- 
plete C/0 nuclear burning in the progenitor WD star. 

In Fig. 19 wc show the Si ii A6355 line profile for 
SN la in the CfA sample for which we have spectra 
prior to — lOd from _B-band maximum light (excluding 
1991T/1999aa-like objects and SN 20Q6cc; see caption) 
. There is a large variation in the maximum Si velocity 
(roughly corresponding to the blue edge of the Si li A6355 

absorption profile), from 28000 km s"^ for SN 2003 W 

to 18000 km s"^ for SN 1998aq. These values are 

consistent with the lower limit of ~ 12000 km s~ found 
by Mazzali et al. (2007) for the outer extent of the Si- 
rich layer in SN la ejecta, although as noted by Branch 
et al. (2007) this layer generally extends to much higher 
velocities, in conflict with pure deflagration models (e.g., 
Maedaet al. 2010b). 

Mazzali et al. (2007) assume that pre-maximum 
Si II A6355 profiles arc affected by high-velocity features 
(see next section) that are not representative of the bulk 
of intermediate-mass elements. In Fig. 19 we clearly see 
variations in the line-profile morphology of the Si II line, 
from single absorptions (both broad as in SN 2003W and 
narrow as in SN 1998aq) to evidence for multiple com- 
ponents (e.g., SN 2006cp, SN20071e, SN 2005cf). In the 
latter case the absorption profile is interpreted as con- 
sisting of a photospheric and a "detached" high- velocity 
component, as in the analysis of SN 2005cf by Wang et al. 
(2009b). In this interpretation both components would 
then be blended together in SN la displaying a single 
broad absorption profile (see, e.g., Tanaka et al. 2008). 

The emission component of the P-cygni profiles in 
Fig. 19 appears blucshiftcd with respect to the rest- 
frame wavelength of the Sin transition (6355 A). In 
SN 2003W, the peak emission occurs at a Doppler veloc- 
ity of ~ —6000 km s~ . Such peak-emission blueshifts are 
routinely observed in Type II supernovae, and result from 
the steep density gradient in the ejecta, confining the ab- 
sorbing and emitting regions (as seen by an external ob- 
server) to a narrow range of impact parameters (Dessart 
& Hillicr 2005). Sizeable peak-emission blueshifts are 
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Fig. 19. — Zoom in on the Si II A6355 line for SN la for which 
WC have spectra — lOd from B-band maximum or earlier. The spec- 
tra are roughly ordered according to the maximum velocity at which 
the Si IIA6355 line absorbs (from ~ -28000 kms"^ for SN 2003W to 
~ —18000 km s^ for SN 1998aq), and arc color-coded according to 
their spectroscopic subclass as defined by Branch et al. (2006). We 
do not show spectra of 1991T/1999aa-like SN which generally have 
weak to non-existent Si II absorptions at these epochs, nor the low-S/N 
highly-reddened spectrum of SN 2006cc. Tentative identification of a 
C II A6580 absorption is indicated for some SN la {downward-pointing 
arrows). The spectrum of SN 1990N (Leibundgut et al. 1991) is not 
part of the CfA sample presented in this paper, but it was obtained 
through the CfA SN Program (sec Table A2). 
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also expected theoretically in SN la (Blondin ct al. 2006), 
although line overlap renders their identification prob- 
lematic (see discussion by Branch et al. 2007). In several 
cases we note a small absorption notch in the Si ii A6355 
emission profile (downward-pointing arrows), usually at- 
tributed to C IIA6580 (see § 6.3), which obviously bi- 
ases the measurement of the velocity at peak emission to 
higher apparent blueshifts. 

6.2. High-velocity features 

The Ca li near-infrared triplet (IRT; 8498 A, A8542 A, 
and A8662A) is less subject to overlap with lines from 
other ions than Ca li A3945 and is hence more appropri- 
ate to search for high- velocity absorption features (HVF) . 
Mazzali et al. (2005a) found HVFs to be present in all 
SN la spectra with early-time data covering the Ca ii 
IRT, with evidence of similar features in the Si ii A6355 
line for some. They suggested that an abundance en- 
hancement in Ca and Si was unlikely to be the sole 
cause, as extreme enhancements were required in the 
outer ejecta (e.g., factor of 20 increase in the Ca abun- 
dance above 20000 km s"^ for SN 1999ee; Mazzah et al. 
2005b). A density enhancement was found to be a more 
likely explanation, either intrinsic to the explosion or 
through an interaction with circumstellar material. In 
the former case one expects a large diversity in HVFs 
due to line-of-sight effects dependent on the geometry of 
the high- velocity material (e.g., thick torus or clumpy 
structure; see Tanaka ct al. 2006). There is direct evi- 
dence for the multi-dimensional nature of HVFs in one 
SN la, SN 2001el, where spectro-polarimetric observa- 
tions revealed the high-velocity component of the Ca ii 
IRT to be polarized at the - 1% level (Wang et al. 2003; 
Kasenet al. 2003). 

Figure 20 shows Ca iiA3945 (black) and IRT (red) line 
profiles in pre-maximum SN la spectra, ordered by phase. 
The diversity in line-profile morphology noted by other 
authors (e.g., Tanaka et al. 2008) is clearly seen, ranging 
from single to multiple detached absorption profiles. The 
correspondence between the Ca ii A3945 and IRT absorp- 
tion profiles is obvious in some cases (e.g., SN 1997bq), 
more questionable in others (e.g., SN 2001da) and nonex- 
istent for a few cases (e.g., SN 2002he, SN 1996X). One 
possible source of confusion comes from the overlap of 
the Si IIA3858 hue with Ca iiA3945, which can lead to 
an apparent HVF in the H&K line where none is visi- 
ble in the IRT (e.g., SN 1996X, SN 1989B). The dotted 
lines in Fig. 20 show the location of maximum absorp- 
tion by Si iiA3858, assuming it has the same absorption 
velocity as Si iiA6355. In several cases it coincides with 
a blue absorption in the Ca ii A3945 absorption profile, 
sometimes resulting in a fiat-bottomed absorption (e.g., 
SN 20021ic). Another difficulty when interpreting HVFs 
associated with Ca ii is that a mass fraction as low as 
^ 10^^ is sufficient to yield a strong Ca ii absorption at 
high velocities (see, e.g., Blondin et al. 2011a). 

Tanaka et al. (2006) noted that if HVFs result from 
a density enhancement (regardless of its origin), other 
lines (e.g.. Si iiA6355) should be affected, and the rela- 
tive strengths and velocities of these HVFs should be cor- 
related in some way. Stanishev et al. (2007) suggest that 
the relative strengths of HVFs in Ca ii A3945, Si ii A6355, 
and Ca ii IRT are indeed correlated in SN 2003du, 



but Tanaka et al. (2008) find that the velocities of the 
Si II A6355 and Ca il IRT high- velocity features in several 
SN la are at odds with this hypothesis. We note based 
on Figs. 19 and 20 that five of the SN la with the largest 
extent of Ca ii IRT absorption (SN 2002bo, SN 2005cf, 
SN 1997bq, SN 2002dj, SN 2003W) also have some of 
the largest Si ii A6355 absorption velocities at early times 
(the same is also true for SN 2009ig, shown in Fig. 20 but 
not in Fig. 19). Interestingly, Tanaka et al. (2006) also 
suggest that a combination of density and abundance en- 
hancements might provide a better explanation for the 
origin of HVFs, on the account that dense blobs covering 
the entire photosphere would result in Si ii A6355 absorp- 
tion velocities in excess of 20000 km s~^, "never observed 
in any SN" . In fact, such high absorption blueshifts are 
seen in one supernova from the CfA sample, SN 2003W 
(see Fig. 15), and more recently in SN 2009ig (Foley et al. 
2012). 

6.3. Evidence for unburnt carbon 

Inferring the presence of unburnt carbon in SN la spec- 
tra can yield important clues on the properties of the 
explosion. Pure deflagration models predict that large 
amounts of unburnt C/0 fuel are left over, sometimes 
mixed to the innermost ejecta regions (see, e.g., Ropke 
2005), although one then predicts strong forbidden lines 
of C I and O l in late-time spectra that are not observed 
(Kozma et al. 2005). Delayed-detonation models, on the 
other hand, consume all but the outermost regions of the 
SN ejecta, leaving varying amounts of unburnt C/0 at 
high expansion velocities (e.g., Maeda et al. 2010b) which 
can affect early-time SN la spectra. 

As noted by several authors (e.g., Tanaka et al. 2008), 
carbon is expected to be once-ionized in the ejecta of 
SN la at early times, such that the "best bet" at con- 
straining the abundance of carbon, let alone inferring its 
presence, is through its strongest optical line C ii A6580 
(although see the analysis by Marion et al. 2006 using C i 
lines in NIR spectra). Until recently, firm detections of 
this line had been reported in only a handful of SN la, the 
most convincing cases being those of SN 2006D (Thomas 
et al. 2007) and SN 2006gz (Hicken ct al. 2007a), both of 
which show gradually fading C ii A6580 lines with respect 
to the earliest detection. The scarcity of C ii detections 
led to the suggestion of a clumpy carbon distribution in 
the outer ejecta of SN la combined with line-of-sight ef- 
fects (e.g.. Branch et al. 2007). 

In a recent survey of C ii features in SN la spectra, 
Parrent et al. (2011) find evidence for C il A6580 absorp- 
tion in ^ 30% of SN la in their sample, a significantly 
larger fraction than previously suspected (e.g., Thomas 
et al. 2007). Moreover, they show that these features 
are most common in SN in the low-velocity gradient 
group, in line with the idea that LVG SN la correspond 
to off-center delaycd-dctonations seen from the deflagra- 
tion side (Maeda et al. 2010a; although sec Blondin ct al. 
2011a for an interpretation of velocity gradients in the 
context of symmetric explosions), or at least to WDs that 
underwent less intense burning (Tanaka et al. 2008). 

Thomas et al. (2011a) investigated the presence of 
C II A6580 in early-time SN la spectra from the SNfac- 
tory collaboration and found definite detections in 5 
SN la. AU five SN are part of the CN/LVG subclass and 
are otherwise unremarkable, apart from SN 2005cf which 
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Fig. 20. — Ca II A3945 (H&K; black) and A8579 (IR triplet, or IRT; red) line profiles in pro-maximum SN la spectra, ordered by phase, and 
plotted in velocity space {bottom x-axis). The spectra have been smoothed using an inverse-variance Gaussian weighted filter with dX/X — 0.003. 
The vertical dotted lines indicate the absorption velocity of the Si II A3858 line {top x-axis) assuming it is equal to that of the Si llA6355 line. 
By comparing the H&K and IRT profiles and taking into account the likely location of maximum absorption by Si II A3858 one can interpret 
absorptions in the blue wing of the Ca II A3945 line as being due to high- velocity ("detached") Ca II or normal- velocity ( "photospheric") Si II, 
although a definitive answer can only come from detailed spectroscopic modeling. Spectra marked with an asterisk are from the CfA SN Program. 
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displays clear high- velocity features in the Si ii A6355 
and Ca ii IRT lines. Of the remaining four SN, three 
(SN 2005el, SN 2005ki, and SNF20080514-002) are char- 
acterized by blue colors and narrow light-curve widths, 
suggesting that carbon-positive SN la may be preferen- 
tially associated with a specific photometric behavior, 
and one (SN 2005di) is highly reddened by host-galaxy 
dust. Thomas et al. (2011a) estimate that ~ 20% of 
"normal" SN la display a C ii A6580 absorption as late 
as — 5 d from _B-band maximum light, in line with the re- 
sults of Parrent et al. (2011) [see also the paper recently 
submitted by Silverman & Filippcnko 2012]. They also 
illustrate using a toy model that a spherically symmetric 
distribution of carbon can account for the C ii A6580 line- 
profile morphology, although they cannot rule out the 
large-scale asymmetries inferred by Maeda et al. (2010a) 
to explain the spectroscopic diversity of SN la. 

In Fig. 21 we show the Si ii A6355 region in previously 
unpublished spectra from the CfA SN Program where we 
tentatively identify an absorption due to C iiA6580. For 
each of the 21 SN la in this figure we only show the ear- 
liest spectrum where the C iiA6580 absorption is seen, 
except for a few cases where the earliest spectra had too 
low S/N (SN 2004at, SN 2004fu, SN 2006ax, SN 2007cq, 
and SN 2008 A). Furthermore, we do not show spectra 
of SN la already included in Fig. 19. In all this repre- 
sents a total of 24 SN la not included in the study by 
Parrent et al. (2011)^3. A C iiA6580 line in SN 2005el 
was already noted by Thomas et al. (2011a) in their ear- 
liest spectrum^^. We do not attempt to confirm these 
detections with parameterized spectrum synthesis codes 
as done by Parrent et al. (2011) using SYNOW (e.g.. 
Branch et al. 2005) or by Thomas et al. (2011a) using 
SYNAPPS (Thomas et al. 2011b), but most of the ten- 
tative C II identifications (17 out of 24) correspond to 
firm detections in spectra from different nights. Of the 
remaining 7 SN la, 3 only have a single prc-maximum 
spectrum (SN 2004bg, SN 2004fz, and SN 2005mz), and 
2 have possible detections in multiple spectra (SN 2000dn 
and SN 2003U). For SN 2001ep, we have a definite detec- 
tion in the — 7d spectrum but none in the —3d spectrum. 
The small notch seen in the —2d spectrum of SN 2007ax 
is no longer visible one day later. 

The majority of SN la with C iiA6580 detections in 
Figs. 19 and 21 are part of the Core Normal subclass 
(19 out of 30). Moreover, 7/9 CN SN la with spectra 
prior to — lOd from _B-band maximum (see Fig. 19) show 
C II absorptions, suggesting that this feature is a generic 
property of this subclass, and highlighting the need for 
very early-time spectroscopy. The other 11 SN la are al- 
most equally divided amongst the SS (4 objects), CL 
(3), and BL (4) subclasses. The carbon-positive SS 
SN la include the 1991T-like SN 2007cq, the 2002cx- 

^^ An even more recent paper by FolatcUi et al. (2012) confirms 
C II detections in 4 SN la from Fig. 21 (SN 2005el, SN 2006ax, 
SN 2007A, and SN 2007af). They argue, however, that SN 2007ax 
shows no such absorption. They also report a C II detection in 
pre-maximum spectra of SN 2008bf. We have spectra of this SN, 
but only past maximum light, by which time the C II absorption 
feature is no longer detected. 

^^ Thomas et al. (2011a) also report a C II detection in 
SN 2005ki, for which we have one FLWO 1.5m-|-FAST spectrum 
taken at +1 d past _B-band maximum light, and SNF20080514-002, 
for which wc have spectra only starting at -|-2d (see Tabic Al), cf. 
— 8d for the latest detection by Thomas et al. (2011a). 



like SN 2008A, and the possibly super-Chandrasekhar 
SN 2006gz. One of the carbon-positive CL (SN 2007ax) 
has an obvious 1991bg-like spectrum, while the remain- 
ing two (SN 2002fb and SN 2005mz; both similar to 
SN 2006D) show evidence for more moderate Ti ii ab- 
sorption. Of the 4 carbon-positive BL SN la, only 
one (SN 2004fu) is part of the high- velocity subclass of 
Wang et al. (2009a), with a high velocity gradient in 
the Si IIA6355 line {v ~ 120 km s^ day" measured be- 
tween -1-2 d and -|-22d). This also confirms the results of 
Parrent et al. (2011) who find only one carbon-positive 
SN la of the HVG group in their sample (SN 2009ig; 
sec also Foley et al. 2012). Parrent et al. (2011) suggest 
line blending as a possible explanation for the apparent 
paucity of high- velocity carbon-positive SN la. 

The C II A6580 blucshifts at maximum absorption in 
this sample of 30 SN la are typically in the range 



-10000 > 



Wabs 



> 



44000 km s" 



with the no- 



table exceptions of the 2002cx-like SN 2008A (wabs ~ 
—8000 km s^ at — 6d, similar to the Si IIA6355 absorp- 
tion velocity at the same phase; see Fig. 15) and the pecu- 
liar SN 2006gz (uabs « -15500 kms"^ between -13d and 
— lid), both part of the SS subclass. The absorption ve- 
locities of C II A6580 and Si ii A6355 are similar for most 
objects at any given pre-maximum epoch, their ratio be- 
ing within 10% of unity (see also Fig. 8 of Parrent et al. 
2011). Two SN la (SN 1990N at -14d and SN 2005cf be- 
tween — 13d and —lid) have Uabs(C ii)/wabs(Si ii) < 0.8, 
but both display high- velocity features in the Si ii A6355 
line (see Fig. 19) which biases this ratio to lower values. 
Conversely, SN 2006gz, SN 2008Z (both SS SN la), and 
the 1991bg-likc SN 2007ax display Wabs ratios closer to 
1.2. 

The maximum velocity at which C ii A6580 absorbs 
(wmax) is in the range —12000 > Wmax f^ — 16000 kms~^ 
for most SN la. Since Wmax and Uabs are strongly 
correlated, we find the same two outliers as above 
(SN 2008A, Umax « -10000 km s"^ SN 2006gz, |w,nax| > 
18000 km s~ ). Within the same sample, the maximum 
velocity at which Si ii A6355 absorbs is typically signifi- 



cantly higher (-20000 



> 



> -28000 km s"'). This 



suggests that some silicon is ejected at larger velocities 
than the bulk of unburnt carbon, incompatible with cur- 
rent delaycd-detonation models (see, e.g., Maeda ct al. 
2010b). As noted by Branch et al. (2007), the pulsating 
reverse detonation models of Bravo & Garcia-Scnz (2006) 
could provide a viable explanation for the low carbon ve- 
locities compared to the bulk of IME. 

Thomas ct al. (2011a) suggest that the presence of C ii 
lines in normal SN la could be preferentially associated 
with blue colors and narrow light-curve widths, based 
on the SALT2 parameters (xi,c). They note, however, 
that this trend does not appear to hold for the sam- 
ple analyzed by Parrent et al. (2011). We have reliable 
SALT2 fits for 24 out of 30 SN la in our carbon-positive 
sample, and do not find evidence for bluer colors with 
respect to the carbon- negative sample (171 SN la with 
rehable SALT2 fits). The median SALT2 xi is approx- 
imately —0.6 for both samples. Since the SALT2 color 
parameter docs not distinguish between intrinsic color 
and extinction by dust, we re-ran this analysis using the 
intrinsic B — V color at B-band maximum and A77ii5(i3) 
values from BayeSN. Within the Core Normal sample. 
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Fig. 21. — Tentative identification of a C II A6580 absorption feature {downward-pointing arrows) in previously unpubUshed prc-maximum 
SN la spectra from the CfA SN Program, ordered by phase. The spectra have been smoothed using an inverse-variance Gaussian weighted 
filter with dX/X = 0.002, and are color-coded according to their spectroscopic subclass as defined by Branch et al. (2006). We do not 
include the spectra with possible C IIA6580 absorption features already shown in Fig. 19. 



the carbon-positive SN la appear to occupy a narrower 
A7TT,i5(i?) range (1 < Ami^{B) < 1.5 mag) than the 
carbon-negative sample (0.8 < Ami5{B) < 1.6 mag), 
but a K-S test shows that the two distributions are con- 
sistent at the p = 0.37 > 0.05 significance level. The 
distribution of intrinsic B — V color is the same for both 
subsamples (p = 0.34). These conclusions hold when 
considering only SN la with persistent C ii features (i.e., 
present later than — 5d), although the numbers become 
then too small for a robust analysis. 



7. NEBULAR SPECTRA 
7.1. Nebular line widths and luminosity 



We have 21 spectra of 12 SN la with spectra taken 
later than 150 d past maximum in the CfA sample pre- 
sented in this paper (see Table 2). At these epochs, 
the spectra are dominated by forbidden lines of Fe ii 
and Fe ill, reflecting the NSE composition of the in- 
ner ejecta regions. Based on the observed correlation 
of the FWHM of the iron emission feature centered at 
~ 4700 A with Ami5{B) in nebular SN la spectra from 
the Asiago SN archive, Mazzali et al. (1998) found a re- 
lation between the kinetic energy of the explosion and 
the mass of synthesized ^^Ni. In particular, the low- 
luminosity SN 1991bg was associated with a low mass 
of ^^Ni but also a low explosion energy, suggesting that 
the fast-evolving light curve was the result of a small 
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Fig. 22. — Relation between decline-rate parameter Ami5(B) 
and FWHM of the Fe emission feature at ~ 4700 A (which in- 
cludes emission from [Fe ll] and [Fe III]) in nebular SN la spec- 
tra. The points are color-coded according to the spectroscopic sub- 
classes defined by Branch et al. (2006). Two SN with no Branch 
class (SN 2006dd and SN 2007sr) are shown as open stars. We 
have used previously unpublished data from the CfA SN Program 
supplemented with our own measurements on publicly-available 
data from the literature, where we have excluded SN 1989B 
and SN 2003cg due to the large associated error (crpwHM ~ 
6000 km s~^ and 3500kms~\ respectively), and SN 2002cx due 
to the low S/N of the spectrum. The dashed line corresponds to 
the original relation published by Mazzali et al. (1998). The inset 
shows a Magellan Clay-|-LDSS3 spectrum of SN 2007sr at 190 d 
past maximum with the Fe A4700 feature highlighted. 

(sub-Chandrasekhar) exploding mass. 

We revisit the relation between the FWHM of the neb- 
ular iron feature at ^ 4700 A and Ami5{B) in Fig. 22 
(see Fig. 2 of Mazzali et al. 1998). We have used our 
own FWHM measurements (see Table 7) on the publicly- 
available data from Mazzali et al. (1998) (i.e., with the 
exception of SN 1989M, SN 1991M, and SN 1993L), sup- 
plemented with 10 of the 12 SN la from the CfA SN Pro- 
gram with nebular spectra (we exclude SN 2002cx due to 
low S/N and SN 2003cg due to the large associated error 
on FWHM) , and additional data from the literature (in- 
cluding a nebular spectrum of SN 1991T from Schmidt 
et al. 1994, publicly available in the CfA SN Archive; see 
Table A2). When several measurements are available for 
a given SN, we plot the error-weighted mean FWHM. 
We also show the original relation published by Mazzali 
et al. (1998) {long-dashed line). 

We have performed multiple least-squares fits of the 
relation Ami^iB) = a + &(FWHM - 15) [FWHM in 
units of lO^kms"^] of Mazzali et al. (1998) to the data, 
the results of which are reported in Table 8, including 
the reduced xt and Pearson correlation coefficient (r). 
We also include the results of fits to the original data 
of Mazzali et al. (1998) (see their Table 1). The en- 
tire CfA-l-Literature sample yields a slope (6 = —0.069) 
that is somewhat shallower than the relation published 
by Mazzali et al. (1998) (& = -0.103), although the cor- 
relation is strong in both cases {r = —0.71 and —0.80, 
respectively). The large xt indicates that the scatter 
about the best-fit linear relation is real and not the re- 
sult of measurement uncertainties. Mazzali et al. (1998) 
suggested that this intrinsic scatter could be taken as 
evidence for a second parameter affecting the light-curve 
widths (in addition to the '^^Ni mass) or the velocities 



of the iron-dominated nebula (in addition to the kinetic 
energy of the explosion). 

As noted by Mazzali et al. (1998), the low-luminosity 
SN 1991bg appears to be a clear outlier, being the 
only SN la with FWHM < 10000 km s"^ (FWHM w 
2500 km s~ ). Excluding this single SN from our sam- 
ple results in a similar slope (b = —0.079) but a sig- 
nificantly degraded correlation (r = —0.37). This is 
somewhat different to the result obtained using the 
original measurements of Mazzali et al. (1998), where 
the exclusion of SN 1991bg leads to a steeper slope 
(6 = —0.169) and a minor impact on the correlation 
strength (r = —0.69). The new CfA-|-Literature sam- 
ple indicates that the scatter about the best-fit linear 
relation (excluding SN 1991bg) is ^ 30% larger than for 
the original Mazzali et al. (1998) sample (~ 50% larger 
when including SN 1991bg). Further excluding the low- 
luminosity SN 1986G (i.e., only considering SN la with 
Ar7ii5(i?) < 1.6 mag), there is essentially no correlation 
between Ami5(B) and the FWHM of the 4700 A Fe line 
(r = —0.17). If the latter quantity is indeed a proxy 
for the kinetic energy of the explosion, as suggested by 
Mazzali et al. (1998), our measurements imply that the 
kinetic energy of the explosion does not strongly depend 
on Ami5{B), and hence on the SN la luminosity, for 
all but the least luminous events. This is expected in 
delayed-detonation models where the WD star is almost 
completely incinerated, as burning of C/0 to Ni or Si- 
group elements releases similar amounts of energy (see, 
e.g., Hoeflich et al. 1995). One will need to see whether 
additional data (in particular for low-luminosity SN la) 
continue to disfavor a strong relation between both quan- 
tities. 

7.2. Nebular line shifts and intrinsic color 

A recent study of nebular-phase SN la spectra by 
Maeda et al. (2010c) revealed that some forbidden 
lines associated with Ni ii and Fe ii displayed up to 
±3000 km s~ velocity shifts about the central wave- 
length, while others did not (this was first noted in mid- 
IR spectra of SN 2003hv by Gcrardy et al. 2007). This 
result was interpreted by Maeda et al. (2010c) in the con- 
text of a simple kinematic model crafted to reproduce 
basic characteristics of an off-center delayed-detonation 
model. These include an inner high-density region of 
stable iron-group elements and ^^Ni, and an outer low- 
density ^^Ni-rich region. The former results from the 
deflagration phase, and hence presents a global velocity 
offset in the case of an off-center ignition, while the latter 
results from the detonation phase, resulting in a spheri- 
cally symmetric distribution of ^^Ni in the outer region. 
Different spectral lines in nebular SN la spectra thus 
trace different zones of the inner ejecta, and the diver- 
sity in velocity shifts for lines tracing the deflagration ash 
(in particular [Fe ii] A7155 and [Ni ii] A7378) then stems 
from viewing-angle effects: the lines are blueshifted when 
the SN is viewed from the ignition side and redshifted 
when it is viewed from the opposite direction. 

Subsequent investigations by Maeda and co-workers 
showed that these nebular line shifts were related to the 
velocity evolution of the Si ii A6355 line around maxi- 
mum light (Maeda et al. 2010a) and to the intrinsic color 
at maximum light (Maeda et al. 2011), suggesting that 



Diversity of SN la Spectra 



27 



TABLE 7 

FWHM MEASUREMENTS OF THE Fe A4700 FEATURE IN NEBULAR SN lA SPECTRA. 



SN 



Branch 
class 



Ami5(B) 
(mag) 



FWHM 
(kms-i) 



Phase(s) 
(d) 



1981B 


BL 


0.91 


0.05) 


13705 ± 790 


+266 


1986G 


CL 


1.65 


0.04) 


12112 ± 521 


+256,+323 


1990N 


CN 


0.95 


0.05) 


15047 ± 1750 


+212,+227,+244,+245,+254,+279,+308,+331 


1991T 


ss 


0.93 


0.05) 


17044 ± 1507 


+259 


1991bg 


CL 


1.94 


0.05) 


2474 ± 317 


+201 


1992A 


BL 


1.41 


0.02) 


14588 ± 519 


+225 


1994D 


CN 


1.37 


0.03) 


15656 ± 717 


+611 


1994ae 


CN 


0.96 


0.04) 


14363 ± 1895 


+152,+366 


1995D 


CN 


1.05 


0.03) 


15612 ± 782 


+283 


1996X 


CN 


1.26 


0.04) 


15039 ± 1459 


+296 


1998aq 


CN 


1.11 


0.04) 


16824 ± 204 


+229, +239 


1998bu 


CN 


1.03 


0.02) 


15473 ± 1577 


+178,+189,+207,+216,+242,+248,+328 


2000CX 


SS 


0.99 


0.02) 


15486 ± 499 


+154,+181 


2001el 


CN 


1.16 


0.03) 


16901 ± 1759 


+397 


2002bo 


BL 


1.08 


0.04) 


20281 ± 35 


+310 


2002dj 


BL 


1.19 


0.05) 


15606 ± 1081 


+219,+271 


2002cr 


CN 


1.23 


0.03) 


13882 ± 720 


+212 


2003du 


CN 


1.07 


0.06) 


15772 ± 1876 


+219,+374 


2003hv 


CN 


1.31 


0.05) 


14433 ± 797 


+322 


2003kf 


CN 


1.10 


0.05) 


18139 ± 1877 


+398 


2004dt 


BL 


1.14 


0.04) 


13294 ± 657 


+347 


2004eo 


CL 


1.38 


0.03) 


14480 ± 630 


+225 


2005cf 


CN 


1.10 


0.03) 


15013 ± 2377 


+317 


2006X 


BL 


1.18 


0.03) 


16653 ± 1754 


+150,+275 


2006dd 


a 


1.07 


0.03) 


15554 ± 1142 


+186,+193 


2007af 


BL 


1.20 


0.02) 


14347 ± 36 


+301 


2007sr 


a 


1.13 


0.06) 


14621 ± 32 


+190 



Note. 
error-weii 



— When several FWHM incasurcnients are available (sec the Phase(s) column), we report the 
fhted mean FWHM. 



^No Branch class for SN 2006dd or SN 2007sr since there are no available spectra around maximum 
light. 
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TABLE 8 
Coefficients of least-squares fits of 
5(S) = a + b(FWHM [lO^kms"!] - 15) TO NEBULAR 
SN Ia spectra. 





Sample a b 


2 

Xi, 


r 


All 

Excl. 

Excl. 


CfA + Literature sample 

1.20 (0.01) -0.069 (0.003) 
91bg 1.21 (0.01) -0.079 (0.004) 
91bg+86G 1.20 (0.01) -0.072 (0.005) 


7.5 
7.4 
6.9 


-0.71 
-0.37 
-0.17 


All 

Excl. 

Excl. 


Mazzali et al. (1998) sample^ 

1.24 (0.03) -0.103 (0.012) 
91bg 1.26 (0.04) -0.169 (0.024) 
91bg+86G 1.25 (0.04) -0.168 (0.034) 


4.4 
2.7 
3.0 


-0.80 
-0.69 
-0.58 



Note. — la errors are given in between parentheses. 
^Measurements taken from Table 1 of Mazzali et al. (1998). 

both latter quantities were also related to viewing angle 
effects in off-center delayed-detonation models. 

We investigate the relation between nebular line shifts 
and intrinsic color in Fig. 23. We use the intrinsic B — V 
color at B-band maximum light derived from BayeSN 
fits as opposed to the intrinsic B"^^^ — V"^^^ pseudo-color 
used by Maeda et al. (2010c, a, 2011) [see also discussion 
in § 5.3]. We were not able to determine reliable in- 
trinsic colors for SN 1986G, SN 2000cx, or SN 2009dc, 
since BayeSN does not include such objects in its training 
sample. We also highlight two objects with an inferred 
visual extinction Ay > 1 rnag, and note that objects at 
redshifts < 0.01 are subject to a larger uncertainty given 
the impact of peculiar velocities on the redshift-based 



distance. We have checked, however, that using external 
distance priors based on Cepheids or Surface Brightness 
Fluctuations results in inferred intrinsic colors that are 
consistent within < la with those reported here. We 
computed nebular line shifts based on a two- Gaussian fit 
to the emission profiles of [Fe ii] A7155 and [Ni ii] A7378. 
The nebular velocity, Wneb, is simply computed as the 
mean velocity of both lines, while the error is set to 
the standard deviation of both measurements (see Ta- 
ble 9). For cases where only one of two lines has a well- 
defined emission profile (usually [Fe ii] A7155), we set the 
error in v^eh to 600kms~ , as done by Maeda et al. 
(2011). Our measurements are consistent within the er- 
rors with those reported by Maeda et al. (2010a, 2011), 
with a scatter about zero difference of '--^ 150kms~^ 
when the same spectrum was used (marked with a "|" 
in Table 9), and of ~ 700kms~ when using a different 
spectrum of the same object. One notable exception is 
SN 2006X for which we infer a ^ 1500 km s~ larger v^ch 
(2832 ± 343kms"^ cf. 1331 ± 164kms"^ reported by 
Maeda et al. 2011). This difference is suspiciously close 
to the host-galaxy redshift (cz = 1550 km s~ from Falco 
et al. 1999) and results from Maeda et al.'s accidental 
use of spectra of SN 2006X available in the SUSPECT^^ 
archive that were already de-redshifted by 1587 km s~ 
(Maeda 2011, priv. comm.; this error does not affect 
their conclusions, however). We include four SN la 
(SN 1994ae, SN 1995D, SN 2003kf, SN 2007af) that are 

^^ http://suspcct.nhn.ou.cdu/~suspcct/ 
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Fig. 23. — Bottom: Intrinsic B ~ V color at B-band maximum 
light, as derived from BayeSN fits to SN la light curves, vs. the neb- 
ular fine shift (tJncb) derived from [Fe ll] A7155 and [Ni ll] A7378 (see 
Table 9), for SN la in the Normal and High- velocity subclasses of 
Wang et al. {2009a). Highly-reddened SN la {Ay > 1 mag) are cir- 
cled in red. Objects part of the low-extinction sample are marked 
with an asterisk (see text for details). We could not determine the 
spectroscopic subclass of SN 2006dd {open star) since no spectra 
are available around maximum light (Stritzinger et al. 2010). We 
have no intrinsic color estimates for the low-luminosity SN 1986G, 
the peculiar SN 2000cx, or the possibly super-Chandrasekhar-mass 
SN 2009dc, but we nonetheless indicate their nebular line shifts 
(crosses). Top: Distribution of tJncb for the Normal [open) and 
High-velocity (filled) subsamples. SN 2006dd and the three SN la 
with no intrinsic color estimates mentioned above are not included 
in either distribution but their names appear in the corresponding 
velocity bin. 

not part of the sample studied by Maeda et al. (2010c, a, 
2011) and for which we have spectra past +150 d from 
i?-band maximum light (see Tabic 2). 

The data points in Fig. 23 are color-coded according 
to their spectroscopic subclass ("Normal" and "High- 
velocity" in the classification scheme of Wang ct al. 
2009a). Our spectroscopic coverage for SN 2007sr only 
starts ~ 2 weeks past i?-band maximum (sec Table Al), 
so we infer its spectroscopic subclass (HV) based on the 
Si II A6355 absorption velocity of — 12300 kms^^ at -f-5d 
reported in CBET 1173 by Naito et al. (2007). We could 
not determine a spectroscopic subclass for SN 2006dd, 
since there are no available spectra around maximum 
light (Stritzinger et al. 2010). 

The distribution of nebular line shifts (Fig. 23, top) 
shows a clear separation at Wncb ~ 1000 km s~ between 
the N/HV spectroscopic subclasses. One notable excep- 
tion is the HV SN 2004dt, whose nebular spectrum was 
shown by Macda ct al. (2010a) to resemble most that of 
the faint SN 1991bg than other HV SN la. We thus repro- 
duce the main result of Maeda et al. (2010a), although 
they base their spectroscopic classification on the velocity 
gradient of the Si li A6355 line defined by Benetti et al. 
(2005) as opposed to its absorption velocity at maximum 
(all the HV SN la in Fig. 23 are classified as HVG by 
Maeda et al. 2010a). In the context of the simple kine- 
matic model of Maeda et al. (2010c), the HV subclass 
consists of SN la viewed from the direction opposite to 



the initial off-center ignition (hence the pronounced red- 
shift of nebular lines), where intermediate-mass elements 
synthesized during the detonation phase are ejected at 
higher velocities. 

The overall correlation between intrinsic B — V color at 
i?-band maximum and nebular line shift is rather mod- 
est (r ~ 0.36; r ~ 0.58 when excluding the peculiar HV 
SN 2004dt). Considering only the 10 SN la part of our 
low-extinction sample (by which we mean SN la with 
E/SO hosts or in the outskirts of spiral galaxies, follow- 
ing Maeda et al. 2011; see SN la marked with an asterisk 
in Fig. 23 and Table 9), the correlation is even lower 
than for the entire sample (r = 0.30). The bulk of SN la 
in our sample with Vneh > —2000 km s^^ (18 out of 21 
with intrinsic color estimates) do nonetheless appear to 
follow a clear trend of redder intrinsic color for larger 
Wncb {r = 0.82). It is unclear whether the two Normal 
SN 1994D and SN 2003hv at -Uncb < -2000 km s~^ sim- 
ply reflect a larger scatter at large nebular line blueshifts, 
or reveal a non-monotonic relation between nebular line 
shifts and intrinsic color (which would then also accom- 
modate SN 2004dt). 

Our results appear to differ slightly from those of 
Maeda et al. (2011), who find a strong correlation and an 
obvious linear relation between intrinsic color and neb- 
ular line shift for their low-extinction sample (which in- 
cludes SN 1994D and SN 2003hv, both modest outliers in 
the relation defined by the bulk of our sample) , and which 
also applies to all other SN la in their sample (including 
SN 2004dt) after excluding three highly-reddened objects 
(SN 1998bu, SN 2002bo, and SN 2006X; their Fig. 4). In 
fact, Macda et al. (2011) use the relation between "intrin- 
sic" ijmax_ymax pgeudo-color^^ and Ami5(i3) found by 
Folatelli ct al. (2010) to correct their own pseudo-color 
measurement. Using this same approach we too find a 
strong correlation between this color residual and v^cb for 
our low-extinction sample (r = 0.72). We note, however, 
that the Ami5(i3)-color relation of Folatelli et al. (2010) 
has a low statistical significance (^ 2cr), and that there is 
no correlation between the ijmax_^max pggudo-color and 
/S.mi^{B) for the low-extinction sample of Maeda et al. 
(2011) (r = 0.05; their Fig. 3). Whether or not we apply 
this correction makes little difference to the strength of 
the correlation between B™'^^ — V'^^'^ pseudo-color and 
I'nob for our low-extinction sample (r = 0.68 with no cor- 
rection). 

In any case these nebular line shifts do seem to provide 
information on the intrinsic properties of SN la indepen- 
dently from Amiz{B). We confirm the lack of correlation 
between Wnob and A77ii5(i?) (|r| < 0.1 whether or not we 
exclude peculiar SN la) noted by Maeda et al. (2010a), al- 
though Maeda et al. (2011) argue that such a correlation 
is expected theoretically given the correlation they find 
between Amir,[B) and viewing angle (for their model 
with 0.3 Mq of 5^Ni; their Fig. 12). The recent study of 
2D delayed-detonation models of Kasen et al. (2009) by 
Blondin et al. (2011a) shows, however, that the relation 
between Ami^^B) and viewing angle is a complex one, 
both non-monotonic and non-linear in character regard- 
less of the ^^Ni mass. While we agree with Maeda et al. 
(2011) that more data are needed to study the relation 



^^ really the apparent B'" 
extinction sample. 



ymax pgcudo-color for a low- 
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between nebular line shifts and AmiQ{B), wc also note 
the limiting predictive power of simple models crafted to 
reproduce the observed trends. 

8. CONCLUSIONS 

We have presented 2603 spectra of 462 nearby SN la 
out to z sa 0.12 obtained through the CfA Supernova 
Program between October 1993 and August 2008 (2065 
spectra of 406 SN la are published for the first time 
in this paper). For comparison, the entire SUSPECT 
archive contains 1741 spectra of 185 supernovae of all 
types, with 867 spectra of 81 SN la (including the 432 
spectra of 32 SN la published by Matheson et al. 2008 
and other spectra from the CfA SN Archive). The SN 
group at Berkeley is also preparing a public release of 
their sample of 1298 spectra of 582 SN la (Silverman 
et al. 2012a), complementary in many ways to the data 
set presented here. The sample of 2447 SN la spectra 
taken with the FAST spectrograph at the FLWO 1.5 m 
telescope constitutes the largest homogeneous supernova 
data set to date. All spectra were reduced in a consistent 
manner (except for 37 MMT spectra from 1993-1997), 
making this a unique sample for detailed spectroscopic 
studies of SN la. The sample spans a large range in 
Ami^{B) (between ^ 0.7 mag and ^ 2 mag) and includes 
spectra as early as — 15d and as late as +611 d from B- 
band maximum. There are 168 SN la with spectra before 
maximum light, and 92 SN la with at least 10 epochs of 
spectroscopy. 

Using this unique data set, complemented with ad- 
ditional data from the literature, we have studied the 
diversity of SN la properties as a function of spectro- 
scopic class, using the classification schemes of Branch 
et al. (2006) and Wang et al. (2009a). We find that 
there is a continuum of properties (i.e., no strict bound- 
aries) between the different subclasses, but that a single 
parameter is unlikely to explain the spectroscopic vari- 
ation between extreme members of each subclass. We 
have studied two hypotheses of Branch et al. (2009), 
namely that (1) SN la from the Core Normal (CN) sub- 
class have a significantly smaller scatter in intrinsic peak 
Mb, and (2) SN la from the Core Normal and Shallow 
Silicon (SS) subclasses age at the same rate that they 
decline photometrically. We reject the first and show 
that the second is not universally true. Last, we exam- 
ine the width-luminosity relation [intrinsic peak Mb vs. 
A77ii5(i?)] for individual spectroscopic subclasses, and 
find that SS— t-CN— i>Broad Line (BL) form a sequence 
of increasingly steeper relations, although a single rela- 
tion is found to adequately describe all subclasses (other 
than the faint CL subclass). The same is true for the 
91T— ^Normal— !>High-vclocity sequence using the Wang 
et al. (2009a) classification scheme. 

We provide updated versions of the relations between 
several spectroscopic indicators (the 7^(Si) ratio of Nu- 
gent et al. 1995 and its pseudo-EW equivalent, the 
pseudo-EW of the Si ii A4130 line, and the 7^(Si, Fe) ra- 
tio of Hachinger et al. 2006) and Ami5(i3). We con- 
firm the strong correlation reported in previous publica- 
tions with a significantly larger sample, but note a sig- 
nificant scatter and possible non-linearities, as well as 
notable outliers, e.g., SN 2006ot in the pEW(Si ii A4130) 
vs. Ami^{B) relation (which we argue is spectroscop- 
ically distinct from SN 2006bt, contrary to claims by 



Stritzinger et al. 2011). 

We then focus on the evolution of the velocity at max- 
imum absorption (wabs) of the defining Si iiA6355 line 
with spectral phase, for the different spectroscopic sub- 
classes of Branch et al. (2006). The CN and BL sub- 
classes display a large variation in Wabs at early times 
(prior to — lOd from _B-band maximum), which is largely 
suppressed at maximum light in the former but persists 
in the latter. The BL subclass has the largest varia- 
tion in the post-maximum evolution of Dabs (the velocity 
gradient). The SS subclass has similar Wabs as the CN 
subclass around maximum light, except for 2002cx-like 
SN la and possibly super-Chandrasekhar events. The 
same applies to the Cool (CL) subclass, with the excep- 
tion of the 1991bg-likc SN 2002es which displays a factor 
of ~ 2 lower Si ii blueshifts (Ganeshalingam et al. 2012). 

We propose alternative measurements of velocity gra- 
dients that do not depend on the phase of the last avail- 
able measurement and can accommodate a non-linear 
evolution of the absorption velocity with time. These arc 
the average decline rate over a fixed phase interval [tp, ii] 
and the instantaneous decline rate at a given phase (i.e., 
in the limit io — >• ii ) ■ The clustering into the three groups 
defined by Benetti et al. (2005) [FAINT, LVG, HVG] 
is less obvious than in their Fig. 3, and the use of the 
instantaneous velocity decline rate at maximum reveals 
a significantly larger range in velocity gradients (from 
- 0kms"May"^ for SN 2005hj to - 400kms"May"^ 
for SN 2003W). 

We revisit the relation of Foley & Kasen (2011) be- 
tween Si II velocity and intrinsic color, re-analyzed by 
Foley et al. (2011) using our velocity measurements. A 
key difference in our approach is the use of intrinsic color 
estimates inferred from the BayeSN statistical model for 
SN la light curves (Mandel et al. 2011), which allows for 
variations in the host-galaxy dust properties on a SN- 
by-SN basis, as opposed to assuming a single reddening 
law for all SN la. The intrinsic colors estimated using 
BayeSN have significantly larger (and probably more re- 
alistic) uncertainties, resulting in a very weak correlation 
between Si ii velocity and intrinsic color (r « —0.20) 
when these are taken into account. SN la from the High- 
velocity subclass have a redder mean intrinsic color that 
those from the Normal subclass, but the difference is only 
at the ^ 2(7 level of statistical significance. Moreover, we 
find that the scatter in intrinsic color is comparable for 
both N and HV samples, contrary to claims made by 
Foley et al. (2011). 

Through a study of early-time spectra, which probe 
the outermost ejecta, we attempt to place constraints on 
the extent and completeness of nuclear burning in SN la. 
Silicon is found to absorb at significantly larger velocities 
(-28000 < Wmax < -18000 kms"^) than the lower limit 
of - 12000 km s"^ found by Mazzah et al. (2007) for the 
outer extent of the Si-rich layer in SN la ejecta, in conflict 
with pure deflagration models. We study the impact of 
high- velocity features (HVF) on the Ca ii near-infrared 
triplet (IRT), and find that SN la with the largest extent 
of Ca absorption correspond to those with the largest 
Si II absorption velocity. By comparing the line-profile 
morphology of the Ca ii IRT and Ca ii H&K absorptions, 
we illustrate the contamination of the latter profile by 
Si IIA3858, which can be confused with an HVF in the 
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TABLE 9 
Intrinsic color vs. nebular line shift. 



SN^ zcMB Wangb Ami5(B) 
class (mag) 



(B - V)o'' 
(mag) 



Phasc<* 

(d) 



V-pc ""Ni 

(kins^-^) (kms"-"-) 



"neb 

(kms^-'-) 



Rcf. 



Notes 



1986G 

1990N* 

1994D* 

1994ae 

1995D* 

1997bp* 

1998aq 

1998bu 

2000CX 

2001el 

2002bo 

2002dj* 

2002er 

2003du* 

2003hv* 

2003kf 

2004dt 

2004eo* 

2005cf* 

2006X 

2006dd* 

2007af 

2007sr 

2009dc 



0.003 
0.004 
0.004 
0.005 
0.008 
0.009 
0.005 
0.004 
0.007 
0.004 
0.005 
0.010 
0.009 
0.007 
0.005 
0.008 
0.019 
0.015 
0.007 
0.006 
0.005 
0.006 
0.007 
0.022 



91bg 

N 
N 
N 
N 

HV 
N 
N 

N 

HV 
HV 

N 

N 

N 

N 
HV 

N 

N 
HV 

N 
HV 



1.65 (0.04) 
0.95 (0.05) 

1.37 (0.03) 
0.96 (0.04) 
1.05 (0.03) 
1.08 (0.06) 
1.11 (0.04) 
1.03 (0.02) 
0.99 (0.02) 
1.16 (0.03) 
1.08 (0.04) 

1.19 (0.05) 
1.23 (0.03) 
1.07 (0.06) 
1.31 (0.05) 
1.10 (0.05) 
1.14 (0.04) 

1.38 (0.03) 
1.10 (0.03) 
1.18 (0.03) 
1.07 (0.03) 

1.20 (0.02) 
1.13 (0.06) 
0.76 (0.01) 



-0.16 (0.06) 
-0.09 
-0.18 
-0.16 
-0.13 
-0.18 
-0.17 
0.06 
-0.15 
-0.05 
-0.07 
-0.13 
-0.21 
-0.14 
-0.12 
-0.07 
-0.12 
-0.15 
-0.09 
-0.21 
-0.13 
-0.09 



(0.04) 


(0.05) 


(0.04) 


(0.06) 


(0.04) 


(0.04) 


(0.05) 


(0.05) 


(0.06) 


(0.05) 


(0.05) 


(0.04) 


(0.06) 


(0.05) 


(0.05) 


(0.04) 


(0.04) 


(0.05) 


(0.04) 


(0.04) 


(0.05) 



+256t 

-1-309 

-l-306t 

-1-152 

-1-283 

-1-300 

-(-239t 

-(-328t 

+146t 

-l-397t 

-1-310 

-l-27lt 

-|-212t 

-|-374t 

-|-32lt 

-1-398 

-1-347 

-l-225t 

-1-317 

-1-275+ 

-l-193t 

-1-301 

-1-190 

-1-281 



-478 

-1204 
-1193 

-1249 
-1567 

834 
2002 
1728 

-644 
-2553 
928 
-2623 
-188 
-123 

2589 

-1696 

583 

1185 
-1190 



-639 



-1065 

-1135 

-172 

1126 

2513 

2313 

582 

-243 

-3097 

1213 

-152 

3075 

-1410 

1079 



-639 
-478 
-2220 
-1204 
-1193 

2539 

-1157 

-1351 

-172 

980 

2258 

2021 

582 

-443 

-2825 

1070 
-2623 
-170 
-123 

2832 

-1553 

831 

1185 
-1190 



±600 

±600 

±220 

±600 

±600 

±410 

±129 

±305 

±600 

±206 

±361 

±414 

±600 

±283 

±384 

±201 

±600 

±25 

±600 

±343 

±201 

±350 

±600 

±600 



This paper 
This paper 

Maedall 
This paper 
This paper 

Maedall 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 
This paper 



Av > 1 mag 
PecuUar SN la 



Ay > 1 mag 
No spectra around max 

Peculiar SN la 



^SN names marked with an asterisk are part of tiic low-extinction sample (see text for details). 
Spectroscopic class of Wang et al. (2009a). SN 2000cx and SN 2009dc arc both peculiar SN la, while SN 2006dd lacks spectra around maximum light. 



^Intrinsic B - 



V color at B-band maximum light from BayeSN. 

"^Rest-frame phase of the spectrum used to determine the nebular line shift, in days relative to B-band maximum. For SN 1994D and SN 1997bp we give 
the phase as reported by Macda ct al. (2010a). A "f" indicates the same spectrum was used as Maeda et al. (2010a, 2011). 

^Mean velocity derived from [Fc II] A7155 (^Fc) and [Ni II] A7378 (i^Ni). The associated error is simply the standard deviation of both measurements. 
When only a single measurement is available, wc set the associated error to 600 km s^ , as done by Macda et al. (2011). We use our own measurements 
when available, and those reported by Maeda et al. (2011) [noted Maedall] otherwise. 



absence of spectra covering the IRT absorption feature. 

We searched for signs of unburnt carbon in prc- 
maximum spectra, a subject that has received much at- 
tention lately with the studies by Parrent et al. (2011) 
and Thomas et al. (2011a), and find evidence for a 
C II A6580 absorption feature in 23 new SN la (4 of which 
have confirmed detections by Folatelli et al. 2012), the 
majority of which are from the CN subclass. More- 
over, 7/9 CN SN la with spectra prior to — lOd show 
C II absorptions, suggesting that this feature is a generic 
property of this subclass, and highlighting the need for 
very early-time spectroscopy. The absorption velocity 
of C II A6580 is in general within ~ 10% of that of 
Si iiA6355, but silicon is present at significantly larger 
velocities than carbon, in confiict with current delayed- 
detonation models. Contrary to Thomas et al. (2011a), 
we find no preference for carbon detections in SN la with 
blue intrinsic colors and narrow light curves, although we 
are currently lacking data to study this association for 
SN la with persistent C ii lines. 

Using nebular spectra of 27 SN la, we investigate the 
relation between Ami5(i?) and the FWHM of the iron 
emission feature at ~ 4700 A found by Mazzali et al. 
(1998), interpreted as a correlation between luminosity 
and explosion kinetic energy, more luminous SN la corre- 
sponding to more energetic explosions. We find the corre- 
lation is largely driven by the low-luminosity SN 1991bg, 
which is clearly separated from the rest of the sample 
in the (Atoi5(_B), FWHM) plane. Further excluding 



SN 1986G from the sample (i.e., only considering SN la 
with Ar7ii5(i?) < 1.6 mag), we find no correlation be- 
tween both quantities (r = —0.17), suggesting that the 
peak luminosity is only weakly dependent on explosion 
kinetic energy for most SN la, as in delayed-detonation 
models in which the WD star is almost completely incin- 
erated. 

Last, we studied the relation between velocity shifts of 
certain nebular emission lines and intrinsic SN la prop- 
erties at maximum light, following the recent work of 
Maeda et al. (2010c, a, 2011), and interpreted in the 
context of an off-center delayed-detonation model. Ex- 
cluding the spectroscopically peculiar SN 2004dt, the 
five SN la in our sample that are part of the High- 
velocity subclass of Wang et al. (2009a) display pro- 
nounced redshifts in the nebular lines of [Fc ii] A7155 and 
[Ni ii] A7378 (wneb > +1000 km s"^) and arc thus con- 
sistent with SN la viewed from the direction opposite 
to an initial off-center ignition. We find a strong cor- 
relation (r — 0.82) of nebular line shifts with intrinsic 
B — V color at i?-band maximum light for SN la with 
I'ncb > — 2000kms~ (18 out of 21 in our sample), in line 
with the results presented by Maeda et al. (2011), but 
note the possibility of a non-monotonic behavior for the 
largest blueshifts. We caution that their toy models do 
not capture the complexity of the variation of ATOi5(i?) 
with viewing angle in off-center delayed detonation mod- 
els, as inferred by Blondin et al. (2011a) based on the 2D 
models of Kasen et al. (2009). 
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The study presented here is far from exhaustive, and 
we have still a lot to learn from the study of SN la spec- 
tra, from constraining the physics of the explosion and 
the nature of their progenitors to determining more pre- 
cise distances on cosmological scales. All 2603 spectra 
will be made publicly available through the CfA Super- 
nova Archive to enable further studies by the rest of the 
SN community. 
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TABLE Al — Continued 



UT Date'' HJD*^ Phase": TcL/Instr."* Range° Disp.* Res.s P.A.'^ |A<J>|' Air.J Flux Std> Sec' Slit™ Exp." Observer(s)° 

W (A) (A/pix) (A) (°) (°) (") (") (s) 

Note. — Tabic Al is published in its entirety in tiic electronic edition of The Astronomical Journal. A portion is shown here for guidance regarding its form and content. 

^ UT at midpoint of obscrvation(s). 

Heliocentric Julian date at midpoint of observation(s). 

^ Rest-frame phase of spectrum in days relative to B-band maximum. For SN la with no reliable estimate for the time of maximum, we indicate the rest-frame days relative to the first spectrum 
preceded by an "@" symbol 

"^ Telescope and instrument used for this spectrum: FAST = FLWO 1.5m+FAST, IMACS = Magellan Baado4-IMACS, LDSS2 = Magellan Clay4-LDSS2, LDSS3 = Magellan Clay+LDSS3, 
MMTblue = MMT-|-Blue Channel, MMTred = MMT-f-Red Channel. 

^ Observed wavelength range of spectrum. 
Spectral dispersion in A per pixel. 

^ Approximate FWHM spectral resolution in A. 
Observed position angle during the observation(s). 

' Absolute difference between the observed position angle and the average parallactic angle over the course of the obscrvation(s). 

J Airmass of the observation. 

^ Standard stars: BD17 = BD-fl7°4708, BD26 = BD+26°2606, BD28 = BD-|-28°4211, BD33 = BD-f33°2642, CD32 = CD-32 9927, EG131 = EG 131, EG274 = EG 274, F15 = Feige 15, 
F25 = Feige 25, F34 = Feige 34, F56 = Feige 56, F66 = Feige 66, F67 = Feige 67, FllO = Feige 110, G191 = G191B2B, H102 = Hiltner 102, H600 = Hiltncr 600, HD19 = HD 192281, 
HD21 = HD 217086, HD84 = HD 84937, HZ44 = HZ 44, HZ14 = HZ 14, L3218 = LTT 3218, L3864 = LTT 3864, L4816 = LTT 4816, vMa2 = van Maanen 2. 



< 
CD 



zn 



™ Spectroscopic slit width. 

" Exposure time. Separate exposures are indicated. 

° Observers: EA = E. Adams, VA = V. Antoniou, HA = H. Arce, JA = J. P. Anderson, ZB = Z. Balog, PBa = P. Barmby, EB = E. Barton, JB = J. Battat, PBc = P. Berlind, GB = G. Bernstein, g. 

WBl = W. P. Blair, SB = S. Blondin, AB = A. E. Bragg, CB = C. Briceno, WBr = W. Brown, NO = N. Caldwell, MC = M. L. Calkins, BC = B. J. Carter, PC = P. Challis, JC = J. R. Cho, ^ 

CC = C. Clemens, ACo = A. Cody, ACr = A. Crook, TC = T. Currie, RC = R. M. Cutri, KD = K. Dendy, AD = A. Diamond-Stanic, ND = N. Dinshaw, JDon = J. L. Donley, JDow = J. J. Downes, O 

KE = K. Eriksen, GE = G. Esquerdo, EF = E. E. Falco, RF = R. Fesen, CF = C. B. Foltz, JF = J. Foster, JGa = J. Gallagher, AG = A. Garg, PC = P. M. Garnavich, IG = I. Ginsburg " 

JGr = J. Graves, NG = N. Grogin, TG = T. Groner, VH = V. Hradecky, HH = H. Hao, CHei = C. Heinke, CHel = C. Heller, JHe = J. Hernandez, MH = M. Hicken, CHi = C. Hill, g 

JHuc — J. P. Huchra, JHug — J. P. Hughes, CHu — C. Hutcheson, RH — R. Hutchins, RJ — R. Jansen, SJ — S. Jha, SKa — S. J. Kannapan, SKe — S. Kenyon, RK — R. P. Kirshner, | , 

DK = D. M. Koranyi, JK = J. Kuraszkiewicz, HL = H. Landt, TL = T. Lappin, NL = N. Lepore, LM = L. Maori, JM = J. A. Mader, AM = A. Mahdavi, EM = E. Mamajek, SMa = S. A. Mao, ^ 

NM = N. Martimbeau, TM = T. Matheson, MM = M. Modjaz, FM = F. Munshi, SMu = S. Muscarella, GN = G. Narayan, PN = P. Nutzman, CP = C. A. Pantoja, BP = B. M. Patten, CC 

KP = K. Ponev, JPe = J. Peters, WP = W. Peters, MP = M. Phelps, JPi = J. Pineda, AR = A. G. Riess, KR = K. Rines, BS = B. P. Schmidt, MS = M. Schrodter, JS = J. D. Silverman, a, 

IS = 1. Song, ST = S. Tokarz, MT = M. Torres, CT = C. Tromonti, AV = A. Vaz, LW = L. Wells, MW = M. Wostovcr, RW = R. J. Wcymann. S. 

P Spectra with ~ 400 A-wide gaps (~6150-6550 A) between blue and red halves. p 

^ The first two standard stars were used to calibrate a spectrum taken with our standard grating tilt, while the third was used to calibrate a spectrum taken with a red tilt. 

^ Spectra accidently ommitted from Matheson et al. (2008). 

^ Spectra strongly affected by dark-current problems following UV flashing, for which we have trimmed off a portion of the spectrum. 
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TABLE A2 
References for spectroscopic data from the literature used in this 

PAPER. 



SN 



References 



1981B Branch et al. 1983 

1984A Barbon et al. 1989 

1986G Phillips et al. 1987, Cristiani et al. 1992=^ 

1989B Wells et al. 1994 

1990N Leibundgut et al. 1991*=, Mazzah et al. 1993=", Gomez & Lopez 1998'' 

1991M Gomez & Lopez 1998=- 

1991T Jeffery et al. 19921^, Schmidt et al. 19941^ 

1991bg Turatto et al. 1996=' 

1992A Kirshner et al. 1993'" , Mazzali et al. 1998 

1994D Patat et al. 1996=' 

1996X Salvo et al. 2001'' 

1997cn Turatto et al. 1998" 

1998bu Gappellaro et al. 2001" 

1999aa Garavini et al. 2004" 

1999ac Garavini et al. 2005", Phillips et al. 2006 

1999aw Strolger et al. 2002 

1999ee Hamuy et al. 2002" 

2000E Valcntini et al. 2003" 

2000CX Li et al. 2001" 

2001ay Krisciunas et al. 2011b 

2001el Wang et al. 2003", Mattila et al. 2005 

2002bo Benetti et al. 2004" 

2002dj Pignata et al. 2008" 

2002er Kotak et al. 2005 

2003cg Elias-Rosa et al. 2006" 

2003du Anupama et al. 2005", Stanishev et al. 2007 

2003fg Howell et al. 2006 

2003hv Leloudas et al. 2009 

2004S Krisciunas et al. 2007" 

2004dt Altavilla et al. 2007 

2004eo Pastorello et al. 2007" 

2005bl Taubenberger et al. 2008" 

2005cf Garavini et al. 2007", Wang et al. 2009b 

2005cg Quimby et al. 2006" 

2005hj Quimby et al. 2007 

2005hk Phillips et al. 2007 

2006D Thomas et al. 2007 

2006X Wang et al. 2008 

2006bt Foley et al. 2010 

2006dd Stritzinger et al. 2010 

2006gz Hicken et al. 2007a'^ 

2006ot Stritzinger et al. 2011 

2007if Scalzo et al. 2010, Yuan et al. 2010 

2009dc Taubenberger et al. 2011 



Note. 



Wo do not include references to SN la spectra from our sample that have previously been published (see references in § 1) 



Spectra dowloaded from the SUSPECT supernova archive. 
Spectra available at the CfA Supernova Archive. 



